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ABSTRACT 

We present DEEP2 galaxy clustering measurements at z ~ 1 as a function of stellar mass, star 
formation rate (SFR), and specific SFR (sSFR). We find a strong positive correlation between 
stellar mass and clustering amplitude on 1-10 h^^ Mpc scales for blue, star-forming galaxies with 
9.5<log(M*/M0)<ll and no dependence for red, quiescent galaxies with 10.5<log(M*/MQ)<11.5. 
Using recently re-calibrated DEEP2 SFRs from restframe _B-band magnitude and optical colors, we 
find that within the blue galaxy population at z ~ 1 the clustering amplitude increases strongly with 
increasing SFR and decreasing sSFR. For red galaxies there is no significant correlation between clus- 
tering amplitude and either SFR or sSFR. Blue galaxies with high SFR or low sSFR are as clustered 
on large scales as red galaxies. We find that the clustering trend observed with SFR can be explained 
mostly, but not entirely, by the correlation between stellar mass and clustering amplitude for blue 
galaxies. We also show that galaxies above the star-forming "main sequence" are less clustered than 
galaxies below the main sequence, at a given stellar mass. These results are not consistent with the 
high-sSFR population being dominated by major mergers. We also measure the clustering amplitude 
on small scales (<0.3 Mpc) and find an enhanced clustering signal relative to the best-fit large- 
scale power law for red galaxies with high stellar mass, blue galaxies with high SFR, and both red 
and blue galaxies with high sSFR. The increased small-scale clustering for galaxies with high sSFRs 
is likely linked to triggered star formation in interacting galaxies. These measurements provide strong 
constraints on galaxy evolution and halo occupation distribution models at z '-^ 1. 

Subject headings: galaxies: evolution — galaxies: active — galaxies: high-redshift 



1. INTRODUCTION 

The non-uniform spatial distribution of galaxies in 
the Universe, known as large-scale structure, likely re- 
flect s the non-uni form spatial distribution of dark mat- 
ter. Kaiser (1984) showed that galaxies are biased trac- 
ers ot the total underlying mass distribution. Measure- 
ments of galaxy clustering therefore reflect, in part, the 
clustering of dark matter at a given cosmic epoch. Fur- 
ther, measurements of galaxy clustering over cosmic time 
constrain the evolving relationship between dark matter, 
which is governed by gravity and reflects cosmological 
parameters, and galaxy properties, which are governed 
by processes associated with baryonic physics. As such, 
galaxy clustering measu rements have the ability to con 



strain both c osmology fPeacock et al. 2001 Seljak & 



Warren 2004| [Eiscnstein et al. 20 05j and galaxy evolu- 



tion physic sl Zhen g et al.|2007||Uo5'roy fe Wechsler|2009 



Zehavi et al.||2U12| ) and are thus a powerful tool m the 
study ot large-scale structure formation. 
The standard method used to measure galaxy cluster- 
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ing is the two-point correlation function (2PCF), which 
quantifles the excess probability of flnding a ga laxy in a 
given volume relative to a random dis tribution (Davis & 



Peebles 
redshitt 



1983 Landy & Szalay 



information, 



1993| . Combined with 
the 2PCt' produces a statistical 
representation of the three-dimensional galaxy density 
distribution as a function of scale. In recent years, the 
clustering of galaxies with respect to dark matter has 
been interpreted using the "halo model" framework. The 
clustering amplitude of both dark matter particles and 
dark matter halos can be analyticall y flt, for a given cos 



mology, using TV- body simulations (Mo & White 1996 
Sheth et"aL]|2001 ). Simple proposed analytic models de- 



tion of scale and mass ( Jing|1998 


IPeacock & Smith'20001 


Seljak||2000 Benson et al.pOOO 


Cooray Sheth 2002 f 



occupation 

galaxies populate dark matter halos using a probability 
distribution generated to match measured galaxy cluster- 
ing statistics, therefore statistically connecting t he bary- 
onic matter of gala xies to the dark matter halos ( Berlind 
&: Weinberg||2002[ ). Such HOD models typically include 
descriptions of the halo structure in terms of the "cen- 
tral" galaxy within a given dark matter halo and the 
sub-halo "satellite" g alaxies that reside wit hin the par- 
ent dark matter halo ([Kravtsov et al.|[2004 ). 

Modern day redshitt surveys are uniquely suited to 
characterize the large-scale galaxy distribution through 
the 2PCF over a wide range of cosmic time, including sur- 



veys of local (z<0.15,|York et al. 2000 


Colless et al.|2001 


Abazajian et al. 


2009), intermediate (z ~ 0.5, Biake et 



al. 


2009 


Eisenstein et al. 12011 iCoil et al.ll20111. 


and 


hie; 


her redshitt Kaiaxies (0.7<z<2, Lilly ct al. 


2007 


van 


Dokkum et al.|2009 


Newman et al.|2012|). Such surveys 



correlat ion exists between stellar mass and halo mass to 
z< 2. iLeauthaud et al. (2012) also used the angular 



with various galaxy properties such as luminosity and 
color, where the clustering amplitude increases for galax- 
ies w ith redder color a nd /or higher l uminosity at least t o 



1 ( |Coil et al.|2"008l hereafter C08; |Zehavi et al.||2011[ ). 
Vhile there is a strong, demonstrable relationsnip be- 



z<] 

While tnere is a strong, aemonstraDie reiationsmp 
tween galaxy luminosity and color with clustering ampli- 
tude, interpreting this relationship in terms of the under- 
lying physics can be complicated, as there are competing 
effects due to varying star formation histories, mass-to- 
light ratios, metallicity, dust, and other complex galaxy 
physics. A cleaner physical relationship in the mapping 
between galaxies and dark matter halos may be derived 
from measuring clustering as a function of properties that 
are easier to interpret, such as stellar mass and star for- 
mation rate (SFR). Using these galaxy properties, one 
may obtain direct constraints on the baryonic processes 
involved in galaxy evolution rather than with intermedi- 
ate dependencies o n lum inosity and color. For exam- 



ple, Zheng et al. (2007) constrained the HOD of the 
luminosity-dependent cliistering results from the SDSS 
and DEEP2 surveys, but the interpretation is compli- 
cated by uncertainties in connecting lumino sity and stel- 



lar mass for different galaxy populations. Zheng et al. 
showed that while one can use the HOD to constrain 
the evolution of the stellar-halo mass relationship, fur- 
ther model refinement requires clustering measurements 
for galaxy samples selected by stellar mass. 

Measuring clustering with respect to galaxy properties 
has been made possible by the advent of larg e redshift 
surveys in the last decade. Using the SDSS, |Li et a\\ 
( [2006^ found that on large scales (> 1 h^^ Mpc), the 
dependence of galaxy clustering on stellar mass closely 
mirrors that of luminosity; the clustering bias is rela- 
tively flat below the characteristic stellar mass scale M* 
and increases exponentially above M* . They further in- 
vestigate how the clustering amplitude depends on other 
galaxy color, morphology, and stellar mass, and find sig- 
nificant scale-dependent differences at low stellar mass 
(log(M*/M0)<lO) between r ed and blue galaxy sam ples. 



At higher redshift {z ^ 0.85), [Meneux et al.| ( |2008[) mea- 
sured the clusterin g of w3200 galaxies in the VVDS (Le 



Fevre et al.| 2003 ) selected by st ellar mass. Similar to 
the local behavior, |Meneux et al.| found that the cluster- 
ing amplitude increases with stellar mass and that the 
clustering amplitude at 9.5<log(M*/M0)<lO.5 is sys- 
tematically lower than in local samples. The clustering 
amplitude at higher stellar masses (log(M*/M0)>lO.5) 
has not evolved in the same span of time within the mea- 
surement errors. 

While it is preferable to measure the 2PCF using spec- 
troscopic redshifts, several studies have used photometric 
redshifts to measure the two-dimensional angular corre- 
lation function with respect to stellar mass. Because pho- 
tometric redshifts depend only on accurate panchromatic 
photometry, larger sample sizes are more easily obtained, 
particularly at higher redshifts where spectroscopy re- 
quires extensive telescope resources. Using data from the 
NEWFIRM Medium Band Sur vey (NMBS, |v an Do kkum 
^' ■ |2011| 



used 

angular clustering measurements to find that a positive 



et al. 2009) between l<z<2. Wake et al 



correlation function aloiig with weak lensing map s from 



COSMOS dCapak et al^|2007l [Ilbert et al |[2009| to fit 

an HOD model between U.2<z<l. i'hey tind that the 
peak of the ratio of M*/Mhaio is a constant with cos- 
mic time, indicating a fundamental connection between 
the fraction of baryons within a given halo mass and the 
conversion rate of those baryons into stars. While an- 
gular clustering though photo-zs is less precise than the 
2PCF, some of the most advanced HOD models to date 
have been built off such data. 

In contrast to stellar mass, there has been very little 
study of ho w galaxy cluster ing depends on SFR. For lo- 
cal galaxies, Li et al. (2008) used SFRs of SDSS galaxies 
to study the clustering dependence on the specific star- 
formation rate (sSFR, the SFR divided by stellar mass), 
though not the dependence on SF R alone. A t highe r red- 
shifts (z ~ 2), a recent study by Lin et al. (2012) mea- 
sured the angular correlation function tor i^zA -selected 
star-forming galaxies in 160 arcmi n^ of the GOODS- 
North field (Giavalisco et al. 2004). Lin et al. found 
a significant large-scale clustering dependence with in- 
creasing SFR and decreasing sSFR. However, no cluster- 
ing analyses using SFR or sSFR have been perf ormed at 
z ~ 1 where the global SFR is rapidly declining ([ H opkins 
fc Beacom||2006l |Zhu, Moustakas fc Blanton|[200^ 



I'he ettects of galaxy clustering on Sl^'R are interest- 
ing for HOD models as the parent dark matter halos of 
galaxies may influence the SFR through various mecha- 
nisms. For galaxy populations that are highly clustered 
and therefore reside in massive halos, the available gas in 
the galaxy may have been stripped, depleting the galaxy 
of the material necessary to sustain further star forma- 
tion. Alternatively, the gas may be heated to the virial 
temperature of the host halo, such that the cooling time 
becomes long enough to suppress star formation. How- 
ever, such star formation "quenching" could arise from 
secular evolution, which would not necessarily be depen- 
dent upon host halo properties. Similarly, the sSFR is a 
measure of the SFR relative to the existing stellar mass in 
a galaxy and reflects how efficiently a galaxy is process- 
ing gas into stars. Determining the clustering properties 
and galaxy bias with respect to SFR and sSFR as a func- 
tion of redshift is therefore crucial to understanding how 
the host dark matter halo and surrounding environment 
of a galaxy play a role in the evolution of star formation. 

While clustering analyses of the 2PCF as a function 
of stellar mass and SFR are relatively rare at intermedi- 
ate redshifts, there are measurements of the relationship 
between these parameters and galaxy environment. En- 
vironment studies generally measure the galaxy overden- 
sity for a given population relative to the mean density at 
the same red shift, often using the N *^ nearest neighbor 
statistic (see Muldrew et al. (2012) for a review of var- 
ious environmental measures). The advantage of using 
the local environment is that it provides better statis- 
tics in what is essentially a cross-correlation between a 
given galaxy sample and the mean density field, rather 
than the auto-correlation that is used in 2PCF analy- 
ses. However, environmental analyses can be difficult 
to model and compare to mock galaxy simulations, in 
part because the mean observed galaxy density is sub- 
ject to selection effects that depend on the particular 



3 



survey. More importantly, density measurements often 
average across a range of physical scales (typically with a 
mean scale of > lh~^ Mpc at z~ 1) to infer correlations. 
Because the "density-defining" population is generally 
measured on scales that are roughly the size of a halo, 
such measurements preclu de smaller sub-halo scales. For 



example, |Li et al.| ([20061 measured the 2PCF of SDSS 
galaxies with respect to their structural parameters and 
found significant clustering dependence on scales below 
<1 h^^ Mpc, a result that was not previously seen in en- 
vironment studies. The scale-dependence of clustering as 
a function of galaxy properties can yield important clues 
as to the physical processes involved in galaxy formation 
and evolution. 

As environment and 2PCF analyses are complemen- 
tary, it is important to compare the trends found in both 
kinds of stu dies to test for consistency, particularly on 
large scales. Cooper et al. (2008, hereafter CP08) used 



data from SUSS and DFEP2 to quantify the evolution 
of the relationships between SFR and sSFR with envi- 
ronment. They found that the relation between sSFR 
and overdensity is essentially unchanged with redshift, 
suggesting that mergers did not drive the decline in the 



global SF R at z<l. Confirming the results from Elbazet 
al] ( |2007[ ), CP08 also found that galaxies with high SFRs 
are tound in high density environments at z ~ 1, which 
leads to an "inversion" of the SFR-density relationship 
between z '-^ 1 and z ^ 0, where high-SFR galaxies are 
found in relatively sparse environments. 

Other studies have measured the evolution of galaxy 
environment with respect to SFR at a fixed stellar mass 
between 0.2<z<l. Several studies report that the av- 
erage color-density relation does not appear to have 
evolved since z<l, indicating that the sSFR-density re- 
lationship is relatively constant with cosmic time since 



z = 1 (Cucciati ct al. 2006; Scodcggio et al.||2009| 


Peng 


et al. 

Tm ' 


2010^ 


Sobral et al.|,2010| |Griitzbauch et al. |2 


Oil). 



CP08, with the caveat that finer sub-sampling of the 
most extreme density regions can reveal a stat istically 



et al 



et al 



significant difference i n galaxy restframe color s (C ooper 
120101). Similarly, [Scodeggio et al.| ( |2009[ ) and [Peng 
pUTDl find that bulk trends in the SFE^density re- 



lation forlield galaxies (typically with log(M*/M0)<ll) 
can also be explained by the ch ange in stellar mass at a 
fixed SFR. Sobral et al. (2010) goes on to suggest that 
most, but iiot all, of the Sl'R-density correlation is due 
to changes in stellar mass in highly star- forming galaxies. 
The general consensus from these studies is that there is a 
fundamental relationship between stellar mass and large- 
scale environment and that correlations between environ- 
ment and other galaxy properties are due to second-order 
effects. 

In this paper, we quantify the relationship between 
dark matter halos and galaxy properties at z '-^ 1 by 
measuring the clustering of DEEP2 galaxy samples se- 
lected by stellar mass, SFR, and sSFR, on scales of 

0.1-20 h^^ Mpc. Our results wiU enable future HOD 
modeling studies at z ~ 1 as a function of these galaxy 
properties and will further be useful for projecting the 
statistical power expected from large-area baryon acous- 
tic osci llation (BAO) survey s at these redshifts (e.g . Big- 



structure of the paper is as follows. In Section 2, we 
present the DEEP2 dataset and derived galaxy proper- 
ties used in this study. In Section |3j we describe the 
galaxy samples defined by stellar mass, SFR, and sSFR. 
Section H] details the clustering analysis methods using 
the 2PCF, while Section[5]presents the results of our clus- 
tering analysis with respect to the various galaxy prop- 
erties. We compare our findings to the recent literature 
and discuss the broader implications of our results in 
Section [6j and we present our conclusions in Section [Tj 
Throughout this work, we assume a ACDM cosmology 
with 17a/=0.3, nA=0.7, (7^=0.8, and h = 0.7. All quoted 
magnitudes are in the AB magnitude system. 

2. DATA 



Th e DEEP2 Galaxy Redshift Survey ( [Newman et al. 
2012 1 is one of the largest, most complete spectroscopic 
redshift surveys of the Universe to z ~ 1. The DEEP2 
survey was conducted on t he Keck II telesc ope using 
the DEIMOS spectrograph ^Faber et al."2003') and mea- 
sured reliable redshifts for 31,656 galaxies from 0<z<1.4 
over four fields dispersed across the Northern Sky. The 
DEEP2 survey targeted galaxies to a magnitude limit of 
i?AB <24.1, and in three of the four fields used a BRI 
color cut to efficiently target galaxies with z>0.7. 

For this study, we initially select galaxies with 
0.74<z<1.4 and a high confidence redshift (>95%, 
zquality >3). As we are interested in measuring clus- 
tering properties as a function of galaxy properties, we 
retain only those galaxies with well-measured values of 
restframe i?-band magnitude (Mb), color {U~B), stellar 
mass (A/*), and star- formation rate (SFR or Because 
the average SFRs are now es timated from the gal axy lu- 
minosity and restframe color ( Mostek et al.|2012a ) , we do 
not exclude red galaxies in our sample, 't hese combined 
selections result in a parent sample of 22,331 galaxies. 

2.1. Calibrations of Stellar Mass and SFR 

The stellar masses used in this study are constru cted 
from the prescrip tion described in_ Lin et al. ( 2007 1 and 
the Appendix of [Weiner et al. 
M7 



{ 20091), whereTE e local 
■ "20031 is 



restframe color-M/L relation (Bell et al 
tended to higher redshift. The prescription depends on 
redshift, restframe _B-band luminosity, and the restframe 
U — B and B — V colors, and is calibrated to a matched 
sample of stellar masses derived from A"-band luminos- 

Figure [l[ 



ity measurements ( Bundy et al 
pares stellar masses trom the D 



2006) 
EF 



com- 

calibration with 



BOSS, Schlegel et al. 2011 Mostek et al. 2012b). The 



those from a matched sample of 7,267 X-band stellar 
masses (K. Bundy, priv. communication) for galaxies 
with 0.74<z<1.4. We confirm that the stellar masses 
used here are in agreement with if-band derived masses 
with a rms scatter of 0.25 dex and a mean offset of -0.05 
dex. The DEEP2 stellar masses are bas ed on a Chabrier 
IMF, in keeping with the [Weiner et al.[ calibration. 

In CP08, SFRs for DEEP2 galaxies were derived from 
measurements of the [O II] emission line and an empirical 
calibration to galaxy i3-band luminosity and SFR per- 
forme d at low redshift (Moustakas, Kennicut & Tremonti 
2006 1 . Applying this local calibration to redshitts z '--^ 1 
poses some unique problems, including accounting for 
luminosity and dust evolution as well as line emission 
from AGN in red sequence galaxies (Yan et al., 2006). 
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Fig. 1. — Comparison of DEEP2 color-M/L stellar masses and 
K-hand stellar masses. The rms scatter between the stellar mass 
estimates is 0.25 dex and the sample mean offset is -0.05 dex for a 
Chabrier IMF. 



More recent work by Mostek et al.l (I2012a| established 
a new empirical SFR calibration for DEEP2 galaxies us- 
ing br oadband optical col or and UV /optic al SED-based 



SFRs (Salim et al. 2009). Mostek et al. showed that 



the strong observed evolution in the resttrame Mb from 
local redshifts 1.3 mag to z = 1) caused the SFRs in 
CP08 to be overestimated by an average of ~ 0.5 dex. 
However, as this evolution causes to first order a zero- 
point offset in the SFR calibration, the relative trends 
between SFR and overd ensity presen ted in CP08 are still 
valid. Here we use the, Mostek et al.l SFR calibration and 



tions commonly used in the literature 


Kennicutt 


1998 


Moustakas, Kennicut & Tremonti||2006 


. 



DEEP2 stellar masses and SFRs described above, using 
log(V'/M*) = log(M*/M0) - \og{tp). We convert the 
stellar masses from a Chabrier IMF to a Salpeter IMF 
by adding a constant 0.2 dex offset and q uote sSFRs 
for a Salpeter IMF. W e note that both the (Moustakas, J 



Kennicut k Tremonti||2006[ ) and [Mostek etl il. (2012a) 
bl^'R calibrations rely on correlated mean trends in large 
galaxy samples, and therefore may not be highly accurate 
on a individual galaxy-by-galaxy basis (typically 0.2-0.3 
dex rms scatter). However, the samples selected from 
such empirical SFR calibrations are statistically repre- 
sentative of the mean SFR in the galaxy population. 

As clustering measurements are more commonly per- 
formed as a function of galaxy luminosity and color, 
we show in Figure [2] how our derived stellar masses, 
SFRs, and sSFRs depend on restframe color and mag- 
nitude, showing contours of each derived property in 
the U — B - Mb plane. The contours of constant stel- 
lar mass generally run perpendicular to the contours of 
constant SFR in this plane. When converted to sSFR, 



these opposing trends cancel out to produce a largely 
luminosity-independent trend with color, reflecting the 
fact that sSFR and restframe color generally reflect the 
light-averaged age of the stellar population. 

3. GALAXY SAMPLES 

In this section, we describe the construction of samples 
for 0.74<z<1.4 DEEP2 galaxies selected by stellar mass, 
SFR, and sSFR. For galaxy clustering measurements, it 
is preferable to use complete samples for which galaxies 
of the same 'type' (selected by magnitude, color, stellar 
mass, etc.) are included to the highest redshift of the 
sample. In this way, within a given sample there is no 
redshift dependence of the galaxy property of interest, 
which would necessitate weighting each galaxy by the 
volume over which it could be observed. Complete galaxy 
samples also facilitate comparison with simulations and 
HOD modeling, where a galaxy population can be well 
described by a few simple selection criteria and a redshift 
limit. 

We create both binned and threshold samples in this 
study. Constructing samples with independent bins fa- 
cilitates identification of clustering trends with respect 
to each galaxy parameter, while constructing threshold 
samples facilitates future HOD modeling applications. 
The binned sample sizes are chosen such that the bin size 
is greater than or equal to the rms scatter in the galaxy 
property and so that there is a sufficient number of galax- 
ies per bin to minimize Poisson errors. For each galaxy 
parameter, we treat the bin with the lowest number den- 
sity, which generally corresponds to the more massive or 
higher SFR galaxies, as a threshold lower limit. 

Our galaxy samples cover two primary redshift ranges: 
0.74<z<1.05 for samples that are complete for both blue 
and red galaxies, and 0.74<z<1.4 for samples that are 
complete only for blue galaxies. The former redshift 
range allows us to compare clustering trends between 
galaxies with different star-formation histories - includ- 
ing both quiescent and star-forming galaxies - while the 
latter allows us to obtain the best statistical measure- 
ments for blue galaxies. 

3.1. Stellar Mass 

As shown in the left panel of Figure [51 galaxy sam- 
ples selected by stellar mass naturally include both red 
and blue galaxies, but the ratio of red to blue galaxies 
strongly depends on stellar mass. In addition, the R- 
band selection of DEEP2 corresponds to a different rest- 
frame wavelength selection with redshift (bluer at higher 
redshift), such that DEEP2 is complete to different stel- 
lar masses for red and blue galaxies as a function of red- 
shift. To probe the full range of stellar masses allowed by 
the data, we construct both color-independent (e.g. 'all' 
galaxies, blue and red) and color-dependent stellar mass 
samples for z<1.05, and we construct m ass samples for 
blue g alaxies only for z<lA. Following Willmer et al. 
( 2006 1 , we use the standard DEEP2 red and blue galaxy 



separation defined by 

{U -B) = -0.032(Ms 



21.62) + 1.035. 



(1) 



Using this color separation, the left panel of Figure |3] 
shows stellar masses for red and blue DEEP2 galax- 
ies as a function of redshift, with boxes indicating our 
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Fig. 2. — The grayscale in all three panels shows the restframe U — B color and B-band magnitude distribution for the 0.74<2<1.4 DEEP2 
galaxy sample, (left) Stellar mass contours span 9.0<log(M,/MQ)<11.5 with levels of Alog(MH,/M0)=O.5. More massive galaxies tend 
to be redder and brighter. The dashed line shows the DEEP2 r estframe gala xv color selection from Equation^ splitting the red sequence 
and blue cloud, (middle) SFR contours calibrated from Mostck et al.J 1 2012a i spanning 0<log(i/>)<1.5 and levels of Alog(i/i)=0.3 Mq yr~^. 
Galaxies with higher SFR are bluer and brighter. SFRs arc correlated m restframe Mb and U — B partially by construction, where optical 
photometry has been used in the calibration, (right) Contours of sSFR spanning — 8.5>log(-(/>/Af»)> — 11 with levels of Alog(i/>/Af*)=0.5 



Galaxies with higher sSFR arc bluer; U — B color and sSFR correlate tightly, as evidenced by the horizontal contours. 
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Fig. 3. — Stellar mass selected galaxy samples, (left) Colored points show all galaxies in the parent sample, where red galaxies are defined 
to be redder than the color cut given in Equation [l| Solid lines show the mass limits for color-independent galaxy samples, where the 
highest mass sample is roughly volume-limited for both red and blue galaxies. Dot-dashed lines show the lower mass sample limits for blue 
galaxies only at a 2<1.05 redshift limit. Red galaxy mass sample limits are given in Tables [T] and [2| (right) Stellar mass bin samples for 
blue galaxies only, defined with finer Alog(M»/MQ)=0.3 bins and approximately complete m DEEP2 to z<lA. The highest mass bin in 
both samples are implemented as a mass threshold in order to allow sufficient statistics in the clustering measurement. 



binned samples to z=1.05. As a starting point, we de- 
fine color-independent samples including both red and 
blue galaxies in two bins in stellar mass. For the high- 
est stellar mass sample, we use a mass threshold of 
log(M*/M0)>lO.8 to ensure a large enough sample size 
such that clustering measurements can be performed 
with reasonable accuracy. This high-stellar-mass sample 
is roughly volume-limited for both red and blue galaxies 
in this stellar mass range. The lower stellar mass color- 
independent mass bin spans lO.4<log(M*/M0)<lO.8 
but is somewhat incomplete for red galaxies at z ~ 1 
due to the R-band magnitude limit in DEEP2. In this 
mass bin, 30% of the galaxies between 0.75<z<0.85 are 
red while only 12% of galaxies between 0.95<2;<1.05 are 
red, indicating that red galaxies are increasingly excluded 
at the highest redshifts. Therefore, we can construct 



only one sample that is complete for both red and blue 
galaxies (with log(M*/M0)>lO.8), and while we quote 
results from the lower mass color-independent bin, we 
refrain from calculating a number density in this sam- 
ple. We further caution that this sample underestimates 
the contribution from low mass, high redshift red galax- 
ies missing in DEEP 2, and we keep this in mind when 
interpreting the results. 

In order to compare the clustering of red and blue 
galaxies separately as a function of stellar mass, we also 
create red and blue mass-selected samples. The mass and 
redshift limits of each sample are given in Table 2, and 
the blue samples are shown in FigurelS] For red galaxies, 
we create two samples with lO.5<log(M*/M0)<ll.O and 
log(M*/M0)>ll.O at z<1.05. To facilitate comparison 
with red galaxies over the same volume, we create three 



6 



blue mass-limited bin samples down to a lower mass of 
log(M*/M0)>9.6, where the blue population is complete 
at z = 1.05. We further create blue mass-limited bin 
samples over a wider redshift range, to z = 1.4, down to a 
mass of log(Af*/M0)=9.9 (right panel of Figure[3|. The 
blue galaxy mass bin samples to z—lA have a slightly 
smaller width (0.3 dex) than the lower redshift samples 
(0.4 dex), but the width is still greater than the rms scat- 
ter in the stellar mass estimates (0.25 dex, see Figure [l]). 
An upper threshold limit of log(Af*/M0)>lO.8 is used 
to define the highest mass bin. 

We also define stellar mass-selected threshold samples 
(not shown in Figure|3| with which to measure clustering 
properties. These threshold samples are useful for com- 
paring with simulations and semi-analytic models, where 
samples may be defined by a given mass limit, and for 
performing HOD modeling. For the color-independent 
stellar mass threshold samples, we only report the num- 
ber density for stellar masses above log(M*/M0)>lO.5 
due to the limited number of red galaxies in the mea- 
sured volume. The sample properties for each stellar 
mass-selected sample are listed in Tables [T] and [2j 



the stellar mass dependence by constructing specific SFR 
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Fig. 4. — SFRs using the |Mostek et al.| ( |2012a l calibration versus 
Mg. Greyscale shows the distribution ol L)l^l^l'2 galaxies between 
0.74<2;<1.4. Contours show lines of constant stellar mass, with 
the same mass levels as in Figure [2] The dashed line shows the 
approximate demarcation between red and blue galaxies, as defined 
in the color-magnitude diagram and Equation nl 



3.2. Star Formation Rate 

As with stellar mass, galaxies at a given SFR can 
have a range of restframe color and magnitude. Fig- 
ure |4] shows the SFRs for DEEP2 galaxies as a function 
of Mb; the dashed line indicates the separation of red 
and blue galaxies following Equation [l] While most red 
galaxies have an estimated SFR below log('0)<O.5 Mq 
yr^^, only those blue galaxies with Mb <-20 have sim- 
ilar SFRs. Figure |4] also shows that SFR samples will 
select across a wide range of stellar mass (shown with 
contours), particularly for blue galaxies where can 
differ by 1.5 dex at a given SFR. Therefore, our analysis 
using SFR-selected samples not only probes the relation- 
ship between clustering and SFR but also has a weak 
dependence on stellar mass. We will explicitly remove 



samples in Section 3.3 

While one could in theory construct color-independent 
samples in SFR given a large enough survey, DEEP2 
does not have a sufficient number of both red and blue 
galaxies at a constant SFR and complete within the same 
volume to obtain a reliable color-independent clustering 
measurement. We therefore create color-dependent SFR 
samples that are approximately complete within their 
restframe color as shown in Figure [5j To mirror the se- 
lection with stellar mass, we choose red galaxies with 
z<1.05 in two SFR samples limited to -0.1<log(V')<0.3 
Mq yr^i and log(?/')>0.3 Mq yr'^. As SFR is strongly 
dependent on M^, particularly in blue galaxies, the 
redshift completeness limits are clear; blue galaxies 
are complete for log(7/')>O.75M0 yr"-'^ below 2:<1.05, 
\og{^)>1.0MQ yr-i for z<1.25, and log(?/')>1.25M0 
yr~^ for z<1.4. The SFR bin sizes are selected to roughly 
reflect the SFR accuracy quoted in Mostek et al. (0.5 dex 
and 0.25 dex rms scatter for red ancTTjlue galaxies, re- 
spectively), while maintaining sufficient sample sizes in 
each bin. The sample properties for each SFR sample 
are listed in Tables [3] and H] 

3.3. Specific Star Formation Rate 

To construct highly-complete samples selected by 
sSFR, we first emphasize that sSFR, as shown in Fig- 
ure [2] is tightly correlated with restframe color. Simi- 
larly, wc find that plotting the sSFR against the rest- 
frame B-band magnitude (Figure pi) shows that nearly 
all DEEP2 blue galaxies with 0.74<2:<1.4 have sSFRs 
between — 10<log(?/'/A/*)>— 8 yr^^, while red galaxies 
have lower sSFRs of log{ip / M^,)<—10 yr~^. Figure [g] 
also shows that like SFR, galaxies at a fixed sSFR have 
a wider range of stellar masses (represented by contour 
lines) in the blue galaxy population than red galaxy pop- 
ulation. Further, a sSFR-selected blue galaxy sample will 
have a residual correlation between the the mean stellar 
mass and luminosity of the sample. This residual corre- 
lation indicates that galaxy luminosity is only a rough 
approximation to the stellar mass within the blue popu- 
lation at z ~ 1. 

Stellar mass and SFR each have different completeness 
limits in DEEP2 as a function of redshift, color, and mag- 
nitude, which could in principle be difficult to disentangle 
when combined as sSFR. However, we can avoid these 
complications by recognizing that sSFR closely traces 
color. We can therefore first limit the parent sample 
by restframe B-band magnitude, Mb-, for red and blue 
galaxies separately, and then define sSFR-selected sub- 
samples free from incompleteness expected at lower lu- 
minosities within each restframe color. The final color- 
dependent sSFR samples will then be complete to the 
magnitude limit of the sample. 

The upper panels of Figure [7] show Mb as a func- 
tion of redshift for red and blue galaxies separately. At 
z<1.05, the parent sample is complete for Mb<— 21 and 
Mb <— 20.5 for red and blue galaxies, respectively. At 
z<1.4 the blue galaxy sample is limited to Mb <— 21. We 
choose not to use a more complicated , color-dependent 
magnitude limit (see Gerke et"ar]|2007 ) so that our sam- 
ples can be easily interpreted in comparison to other 
sSFR clustering studies. 
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Fig. 5. — Galaxy samples defined by star formation rate, (left) SFR versus redshift for all red galaxies in the parent sample. SFR-selected 
bin samples for red galaxies are shown with dashed lines; these samples are complete at z=1.05 for log(i/>)>0.3 and approximately complete 
for log(i/))>-0.1. (right) SFR versus redshift for all blue galaxies in the parent sample. SFR-selected bin samples for blue galaxies complete 
to 2=1.05 are shown with red, dashed lines, while samples complete to either 2=1.25 or 2=1.4 are shown with black, dot-dashed lines. As 
with stellar mass, the highest SFR samples are treated as thresholds. 

each of the galaxy samples. In a given volume, ^(r) mea- 
sures the excess probability of finding galaxy pairs as a 
function of separation over that of a random distribu- 
tion. The excess proba bility can be calculated using the 
Landy & Szalayl ( 19931) estimator, 
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Fig. 6. — Galaxy specific SFR versus the restframe S-band mag- 
nitude, Mg. Contours correspond to constant stellar mass, with 
the same levels shown as in Figure [2] 

After making these initial cuts on Mg and redshift, 
we can construct sSFR samples that are complete to the 
given Mb limit. The bottom row of Figure [7] shows our 
sSFR-selected samples. In general, the lowest luminosity 
galaxies have the highest sSFRs. We define upper limits 
on sSFR of log(V'/M*)<-9.9 M© yr"^ for red galaxies 
and log('0/Af*)<— 8.6 M© yr~^ for blue galaxies. We also 
limit the blue galaxy population at z<1.05 and z<lA 
in order to facilitate comparisons with the stellar mass 
and SFR-selected samples. The individual sSFR sample 
properties are given in Tables [5] and |6] 

4. ANALYSIS METHOD 

4.1. The two-point correlation function 

The clustering analysis performed in this study mir- 
rors that of COS, and we summarize the methodology 
here. We measure the clustering of DEEP2 galaxies by 
calculating the two-point correlation function, ^(r), in 
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Here RR, DD, and DR are the random-random, data- 
data, and data-random pairs, while ur and are the 
mean number densities of random points and galaxies 
in each sample. Each pair is separated by a three- 
dimensional distance r. As such, £^{r) reflects the three- 
dimensional clustering, which must be inferred through 
measurements of the angular separation between galax- 
ies on the sky and a distance measurement along the line 
of sight. 

The redshift of a galaxy reflects both its distance (as- 
suming a given cosmology) and its peculiar velocity. Due 
to the latter effect, the observed clustering signal has 
additional power that appears in redshift space but is 
not attributable to clustering in real space. On small 
scales, the peculiar velocities of individual galaxies within 
collapsed overdensities results in an "elongation" of the 
clustering signal along the line of sight, known as the 
"Fingers of God". On larger scales, the coherent in- 
fall of galaxies streaming into potential wells causes a 
"squashing" of the signal along the line of sight. To re- 
move these redshift-space distortions from the real-space 
clustering signal, we separate the correlation function us- 
ing two separate coordinate vectors. We define two vec- 
tors, s—{vi~V2) for the redshift-space separation and 
Z=^(i;i+i;2) for the mean coordinate of the pair, and 
calculate 

^ = 4^^ (3) 
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Fig. 7. — Galaxy samples defined by specific star formation rate. (Top row) Restframe i3-band magnitude for red (left) and blue (right) 
galaxies in the parent sample are shown as a function of redshift. To 2=1.05, the red galaxies are complete to Mfl=— 21 while blue galaxies 
are complete to Mg=— 20.5. To 2=1.4, blue galaxies are limited to Mg=—21.0. (Bottom row) From the magnitude-limited samples defined 
in the top row, we create sSFR-selected samples with log(i/i/M»)<— 9.9 for red galaxies and log(i/i/A/»)<— 8.6 for blue galaxies. 



We then can express the three-dimensional ^(r) as a two- 
dimensional projection, ^(rp, tt), where rp is the distance 
across the line of sight and tt is the distance along the 
line of sight. We apply the estimator in Equation [2] to 
^(rj3,7r) by counting galaxy pairs. The random samples 
are drawn from the measured selection function for each 



galaxy sample spatial distribution (see Section 4.2 1. 

In order to recover the correlation function signal 
without redshift space distortions, we project the two- 
dimensional correlation function ^(?'p,7r) onto the rp co- 
ordinate across the line of sight. The projection requires 
integrating along the line of sight, following |Davis fc Pee-| 
bles (119831, 



UJp{rp) 



d7r^(rp,7r) 



dy^ (r 



,2\l/2 



(5) 



where y is the real-space separation along the line of 
sight. By modeling ^(r) as a power law, S,^{r /ro)^'' , 
Equation [5] has an analytical solution, 



roV r(l/2)r[(7~l)/2] 

r(7/2) 



(6) 



which can be evaluated for given values of the clustering 
length. To, and correlation function slope, 7. Follow- 
ing COS, we integrate the projected correlation function 
ujp{rp) to a limit of TTmax = 20 Mpc, beyond which 



the DEEP2 sample is not large enough to measure an 
accurate signal. Implementingthe TTmax limit requires 
that we also correct Equation^ to be truncated to the 
same integration limit. However, the resulting ujp{rp) is 
accurate for measured scales where Tp/ir^^y^ < 0.25 (see 
COS for details). 

We fit for a power law of the form ^=(r'/ro) ^ on scales 
of 1-10 ft-^^ Mpc. As we will show, these scales mitigate 
the effects of non-linear behavior at smaller Vp and avoid 
noisy data at larger scales. The final values of tq and 7 
are obtained by performing a minimization between 
the measured ujp{rp) and Equation [6j after correcting for 
the truncation of TTmax- 

For each sampled scale, galaxy pair counts are aver- 
aged over 10 separate DEEP2 pointings and errors on 
ijjp{rp) are computed using the standard error among 
these pointings. While formal errors for tq and 7 are 
computed from the fit, we report the rms errors for 
each parameter using 10 jackknife samples of the data. 
The jackknife error estimates on rp and 7 encapsulate 
the cosmic variance inherent within the DEEP2 data and 
take into account the covariance among the rp bins. Be- 
cause some DEEP2 pointings are correlated with each 
other on the sky, the cosmic variance may be slightly 
unde restimated. Ba sed on calculations with QUICKCV 
( (Newman fc Davis]|2002 ), the actual rms error due to 
cosmic variance could be up to ^ 6% larger than the es- 
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timates from the 10 jackknife samples. We do not include 
this estimated error contribution in our analysis, elect- 
ing instead to work with the error given directly from the 
data. 

In order to facilitate interpretation of the results and 
to avoid relying solely on power-law fits which may not 
describe the data well over all scales, we further esti- 
mate the galaxy bias for each sample. The bias is a 
measure of the clustering of galaxies relative to that 
of dark matter at the same redshift. We generate a 
dark matter ^(rp, n) for t he mean redsh i ft of e ach galaxy 
sample using code from Smith et al. (20031, integrat- 
ing to the same TTmax as in the data. We then calcu- 
late the average bias over scales of 1-10 Mpc us- 
ing 6^=[wj,(rp)]gai/[a;p(r-p)]dark matter and weighting each 
scale equally in the average bias. The average scale of 
the mean bias is fp—A.l h^^ Mpc. The dark matter 
correlation function assumes a ACDM cosmology with 
^^m(0)=0.3 and CT8=0.8. A different assumed value of 
tTg will linearly scale the bias value (e.g. a CT8=0-9 will 
increase the absolute bias by ^ 12%). The bias error 
is estimated from the mean bias rms measured from 10 
jackknife samples, similar to the errors on tq and 7 using 
power law fits. 

4.2. DEEP2 Masks and i/V^ax 

When calculating ^(r^, tt) for a given galaxy sample, we 
generate random catalogs of uniformly distributed points 
spanning the same volume as the data. There are sev- 
eral observational selection effects that must be taken 
into account in the random catalog in order to match 
the spatial distribution of the data. Examples of these 
effects include the overall survey geometry, bright star 
masks, incompleteness due to DEIMOS slit design and 
placement, and variations in the redshift success rate due 
to observing conditions for each slitmask. The effects 
of survey geometry, data gaps, and the varying redshift 
success rates are all accounted for in the DEEP2 window 
function. 

To remove the window function and any redshift se- 
lection effects, we generate a selection function from the 
galaxy redshift distribution and their estimated comov- 
ing distance in a ACDM cosmology. We smooth and 
spline fit the one-dimension spatial distribution such that 
local overdensities due to cosmic variance are removed 
while preserving the overall shape of the distribution. 
The galaxy distributions are calculated in comoving bins 
of Ar =0.011 Mpc, and smoothing of the spline fit 
occurs over scales of Ar=0.03-0.05 Mpc. All selec- 
tion functions are checked by eye such that no power on 
smaller scales is retained in the smoothed selection func- 
tion, which could reduce inferred correlations on those 
scales. Random samples are then generated from the 
probability distribution of the smoothed selection func- 
tion. The procedure ensures that the computed corre- 
lation function statistics are not affected by the survey 
design or completeness of our selected galaxy samples 

Although the correlation function statistics are robust, 
a systematic error is still imprinted on small scales where 
the correlation function is more likely to be incom plete 
due to slit collisions. Using DEEP2 mock catalogs (Yan 
et al.||20d4 |, COS found that the bias varies smoothly as 
a function of scale and varies in size from 25% on the 
smallest scales to 2% on the largest scales. We apply a 



correction for this systematic bias directly to the mea- 
sured ^(rp,7r) and Wp(rp), and we include an additional 
rms error term for the correction which is measured from 
the mocks and added in quadrature to Up{rp). 

To aid in the application of our results to HOD models, 
we calculate the number density for each of our com plete 
sarnples using the nonparametric 1/V^iax met hod ( Fel- 
Itenll 976 ) employed in previous DEEP 2 studies (Blanton 



2006 ; Zhu, Moustakas & Blanton 2009 ) . The volume over 



which a given galaxy can be observed is given by 



- / dVL 



dz 



(7) 



where Dr jz) is the co moving distance and a spatially flat 
universe ( Hogg|1999 ), dfl is the solid angle of sky covered 
in the survey, and f{z) is the probability that a given 
galaxy is targeted and a successful redshift is produced 
in a given DEEP2 mask. We calculate dfl using the same 
DEEP2 masks as those used in our clustering analysis, 
which excludes regions of bright stars and accounts for 
mask edges, and we find that the integrated survey area 
covered in our slit masks is 2.54 deg^. 

The DEEP2 redshift success rate, f{z), is a function 
of both i?-band magnit ude and restframe co lor; these 
are discussed i n detail in Newman et al.| ( 2012 1 (see also 



Blanton|[2006l IWillmer et al.pOO(jl and [Zhu, Moustakas] 
fc Blanton||2009 ). For each observec 



served galaxy, we generate 
1200 simulated galaxies with similar observational prop- 
erties between Zmin and Zmax and calculate fc-corrected 
magnitudes for each simulated redshift. We then apply 
the selection function /(z) to these simulated galaxies 
and produce a fraction which is multiplied by the co- 
moving volume to produce Knax- 

Once Vinax is obtained, we calculate the effective num- 
ber of galaxies in the comoving volume. 
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and divide by the comoving volume between Zmin and 
•Zmax to derive the number density, n. Upper and lower 
limits for n are computed by assuming a Poisson error 
distribution with approximations provided by [Gehrels] 
(1986) . However, because DEEP2 is a small area sur- 
vey, cosmic variance (ctcv) will be a sig nificant source of 
error in the numb er density. Following Zhu, Moustakas] 
& Blanton (2009), we calculate the variance in number 
density trorn each of the four DEEP2 fields and report 
the estimated Ccv along with each of our samples. In 
selected stellar mass cases where DEEP2 is complete, we 
have cross-checked our number densities with the pub- 
lis hed galaxy mass fun ction between 0.75<2<1.0 given 



Bundy et al. (2006) and find good statistical agree- 



ment with their measurements 

5. GALAXY CLUSTERING RESULTS 

5.1. Clustering Dependence on Stellar Mass 

Results for the two-dimensional correlation function, 
^(rp,7r), as a function of stellar mass are shown in Fig- 
ure[8]for the 0.74<z<1.4 blue galaxy samples. The figure 
shows contours of constant correlation strength for three 
bins in stellar mass ranging from 9.9<log(M*/M0)<lO.8 
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Fig. 8. — Contours of constant correlation strength for the two-dimensional correlation function, ^{vp, it), for blue galaxies with 0.74<z<1.4 
binned by stellar mass. The contour lines correspond to correlation amplitudes of 0.25, 0.5, 0.75, 1.0, 2.0, and 5.0, with ^(rp,7r)=l shown 
as a thick line. The number of galaxies in each hin is as follows: 4774 (upper left), 4279 (upper right), 2435 (lower left), and 1134 (lower 
right). The contours have been smoothed by a 1 X 1 Mpc boxcar for clarity in the figure. 



and a fourth sample with a lower stellar mass threshold 
of \og{M^/MQ)il0.8 (lower right panel). In the lowest 
stellar mass bin, (,{rp,7r) is relatively symmetric both 
across and along the line of sight, while at higher stellar 
masses there is an increasing asymmetry in ^(rp,7r). At 
all masses there is evidence of a compression along the 
line of sight on large scales, due to coherent infall into 
gravitational potential wells known as the "Kaiser effect" 
(Kaiser|1987| . Redshift -space distortions on small scales 
( "1" 'ingers ot God" ) also increase dramatically with stel- 
lar mass, indicating that high-stellar-mass blue galaxies 
have greater virial motions and therefore are in systems 
with larger halo mass. These results mirror those of COS, 
who found a similar increase in redshift-space distortions 
with increasing luminosity (see their Figure 7). 

To mitigate the effects of these redshift-space distor- 
tions in our clustering analysis, we integrate the two- 
dimensional correlation function along the line of sight 
to project the clustering signal onto the Vp axis. This re- 
sults in the projected correlation function, a;p(rp), which 
we show for color-independent mass threshold samples in 
Figure 9|a nd color-dependent mass threshold samples in 
Figure 10 We also plot the proje cted dark matter corre- 
lation runction ( Smith et al.|2003[ ) generated for a ACDM 
universe ((78=0.8) at a mean redshift similar to that of 
the galaxy samples (z = 0.9). In all cases, ujp{rp) is 



consistent with a power law on scales larger than >1 
h"^ Mpc where the clustering of galaxies in independent 
halos (e.g the two-halo term) dominates the correlation 
function. 

On scales of rp<0.3 h^^ Mpc, we find a departure 
from a single power law for the highest stellar masses in 
the color-independent sample (Figure [9]) and at similar 
masses in red galaxies (left side of Figure 10 1. The de- 
parture on small scales is seen more clearly m Figure [TT] 
where we divide out the large-scale power law fit over 
1-10 Mpc and plot the relative correlation function 
for red and blue galaxy mass samples. For clarity, the 
figure shows these relative correlation functions with an 
offset of 0.3 dex between different stellar mass samples. 

To determine the significance with which the observed 
clustering on small scales is a departure from the large- 
scale behavior, we perform a t-test between the rela- 
tive large-scale and small-scale clustering amplitude (Fig- 
ure ITT]). For each jackknife subsample of a given galaxy 
sample, we fit a power law over large scales with rp=l- 
10 Mpc and ratio the correlation function averaged 
over rp<0.3 Mpc to the best-fit power law averaged 
over the corresponding small scales. We then compute 
the mean and standard deviation among the 10 jack- 
knife sample ratios and compute a t score in which the 
null hypothesis ratio value is 1. Finally, we compute the 
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Fig. 9. — The projected correlation function, ajp(rp), for color- 
independent galaxy samples selected by stellar mass thresholds and 
limited to z<1.05. The plot shows the two highest mass threshold 
samples with log(Af*/MQ)>10.5; the highest mass sample is com- 
plete for both red and blue galaxies and the lower mass sample is 
slightly incomplete for red galaxies beyond 2:^0.9. The projected 
clustering of dark matter (so lid, red line) is generated from the 
prescription of |Smith et al.|20 03 and shown for a redshift of 2=0.9, 
which is similar to the mean redshift of the galaxy samples. 



corresponding p- value drawn from a ^-distribution with 9 
degrees of freedom. Assuming that there is minimal cor- 
related error between rp<0.3 Mpc and > 1 
Mpc, this procedure encapsulates the error associated 
with both the large-scale power law fit and the covari- 
ance of LOp{rp) on small scales. 

For red galaxies with lO.5<log(M*/Af0)<ll.O, the t- 
test produces a p- value of p=0.005, meaning that a small- 
scale deviation as large as that observed or larger would 
occur only 0.5% of the time if there were in fact no true 
deviation. Red galaxies with log(M*/MQ)>11.0 have a 
small-scale deviation from the large-scale power law fit 
withp = 0.042, a significant deviation at thep < 0.05 sig- 
nificance level. The lower significance in the higher mass 
galaxy sample is due in part to larger correlated varia- 
tion in the large-scale power law slope 7 in the jackknifed 
samples. For blue galaxies (right side of Figure 11), we 
find no statistically significant increase in clustering on 
small scales relative to large scales at any stellar mass 
probed in this study. The difference in the relative cor- 
relation function amplitude between red and blue galaxy 
samples on small scales could indicate a change in the 
central-to-satellite fraction and/or a change in the halo 
mass distribution in the halo mass function ( Zheng et al 



2009 1 . We will further discuss the physical interpretation 
of this result in Section WM 

The left panel of Figure [12] shows the measured corre- 
lation length, ro, when fitting u)p(rp) with a power law 
model of £, = {r/rQ)~^. As the slope of the correlation 
function, 7, can be covariant with ro, we fix 7=1.6 to 
study the trends in rg with stellar mass. We find a strong 
positive correlation between tq and stellar mass for 
log(M*/M0)<ll.O for blue, red, and color-independent 
galaxy samples. At ll.O<log(M*/M0)<11.5, the clus- 



tering length appears to level off at rg ~ 5 Mpc for 
all colors and does not depend on stellar mass, within the 
errors. The tq results could indicate that high-stellar- 
mass blue galaxies at z>l are precursors to similar mass 
red galaxies at z<l, as they likely reside in similar mass 
halos. We note that the DEEP 2 sample has poor con- 
straints for all stellar masses above log(M*/M0)>11.5; 
the survey volume is not large enough to robustly sample 
the rarest, most massive galaxies. The measured values 
of To, 7, and the absolute bias for each stellar mass sam- 
ple are given in Tables [T] and [2] By design, the absolute 
bias has a similar trend to tq where the bias increases 
with stellar mass. We find no significant trend between 
the best fit correlation slope 7 and stellar mass. 

For the threshold samples in Table [T] we also esti- 
mate the minimum and mean halo mass using the mea- 
sured large-scale galaxy bias and the halo mass func 



1999 Jenk- 



tion for centra l galaxies ( Sheth fc Tormen^ 

ins et al.||2001 ), modified to accommodate tne increased 
frequency of satellites a t low halo mass (see the Ap- 
pendix of Zheng et al. (2007) for details). For each 
galaxy sample, we use the halo mass function calcu- 
lated for the mean redshift of the sample in a ACDM 
cosmology ((78=0.8), and we assume a satellite f raction 
of ~ 17% as measured for DEEP2 L* galaxies (Zheng 



et al. 2007). We caution that the minimum and mean 
halo masses are intended as estimates in lieu of a full 
HOD model fit to u)p{rp), which we reserve to a future 
study. Table [l] shows that the mean halo mass ranges 
from <log(A/haio)>= 12-3 h~^MQ in the lowest stellar 
mass blue galaxy sample to <log(Mhaio)>'~^ 13 h~^MQ 
at the highest stellar masses. This mass range is par- 
ticularly interesting as halo mass assembly models pre- 
dict that galaxies at z—1 transition from slow cluster- 
ing growth below 1O^^M to near exponential clustering 



growth abo ve 10^^ Mq (Moster et al. 



2010 



Leauthaud 



et al. 2012). Although the halo masses presented here 



are estimates, we expect that they will help constrain 
the stellar-mass/halo-mass relationship in this important 
mass range at z ~ 1. 

5.2. Clustering Dependence on Star Formation Rate 

We now investigate how clustering properties depend 
on SFR u sing the SFR-selected galaxy samples described 
in Section [3?2| Because complete samples of red and blue 
galaxies probe vastly different SFR ranges, we consider 
only galaxy samples separated by color (see Figure [5]). 
Following the same procedure as with stellar mass, we 
investigate trends in the redshift space distortions for 
blue galaxies as a function of SFR by studying ^(rp, tt). 
We find that blue galaxies at all SFRs display "Kaiser 
infall" on large scales, and all samples have relatively 
small "Fingers of God" on small scales. This implies that 
large-scale infall dominates the ^(rp, tt) signal in our SFR 
samples and contrasts with the significant redshift-space 
distortions seen as a function of stellar mass. The lack of 
"Fingers of God" indicates that the SFR samples must 
span a broader range of stellar masses than mass-selected 
samples, thus diluting the strong "Fingers of God" seen 
in the latter. 

Turning to the projected correlation function, Fig- 
ure |13 shows Wp(rp) for two SFR threshold levels rang- 



ing from -0.1<log(V')<0.4 Mq yr ^ for red galaxies and 
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Fig. 10. — The projected correlation function, ujp{rp), for color-dependent galaxy samples selected by stellar mass. Stellar mass threshold 
samples for red galaxies at 2:<1.05 are shown on the left, while stellar mass threshold samples for blue galaxies at z<lA are shown on the 
right. 
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Fig. 11. — The relative deviation in ujp(rp) from the best-fit power law for galaxy samples separated by red (left) and blue (right) 
restframe color and binned by stellar mass. The power law is fit to ajp(rp) on scales of 1-10 Mpc. Each sample is offset by an arbitrary 
0.3 dex for clarity in the figure. The mean deviation between the relative small-scale clustering below rp<0.3 Mpc and large-scale 
clustering between 0.4<rp<10.5 Mpc is significant in red galaxies at the p < 0.1 level for 10.5<log(M«/MQ)<11.0 and p < 0.001 level 
for log(Af»/M0)>ll.O. We find no significant increase in small-scale clustering for any of the blue galaxy mass samples beyond the best-fit 
large-scale power law. 



three SFR threshold levels between O.75<log(-0)<l-25 
Mq yr^^ for blue galaxies. The samples shown are re- 
stricted to z<1.05 for a direct comparison between red 
and blue galaxy samples within the same volume. We 
find that the clustering amplitude changes much less as 
a function of SFR within the red or blue sample com- 
pared to the difference in amplitude between the red and 
blue samples. Blue galaxies with higher SFR have higher 
clustering amplitude, while there is no detected difference 
in clustering for the two red galax y sa mples as a function 
of SFR, within the errors. Figure [T4| shows the deviation 
of ujp{rp) from a power law fit on scales of 1-10 Mpc 
for red (left) and blue (right) galaxies binned by SFR. 
For all samples with log(i/')<1.5 Mq yr~^, there is no 



significant deviation from a power law on small scales 
(rp<0.3 h-^ Mpc). For the highest SFR blue galaxies 
with \og{ip)>1.5 Mq yr~^, however, we detect a factor 
of seven increase in the mean clustering amplitude below 
rp<0.3 Mpc relative to the large-scale power law fit. 
Performing the same ^-test as was done for stellar mass, 
the excess clustering signal has a p-value of p = 0.038 
and therefore is significant at the p < 0.05 level. 

The measured tq values for color-separated galaxy 
samples binned by SFR are shown in the center panel 
of Figure [12] The clustering scale length increases with 
increasingsFR for blue galaxies and is constant for red 
galaxies as a function of SFR, within the measurement 
error. However, blue galaxies with the highest SFR have 
clustering amplitudes similar to red galaxies, suggest- 
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Fig. 12. — The clustering scale length, rg, measured by fitting a power law to Ljp{rp) over linear scales of 1-10 Mpc. For this figure, 
the power law slope is fixed to 7=1.6 so as to remove degeneracies laetween rn and y . The scale lengths are measured for galaxy samples in 
stellar mass (left), SFR (center), and sSFR (right) as described in Sections |3. 1)3. 3| Vertical error bars are calculated from the rms scatter 
of ro fit from 10 jackknifed sub-samples in each data bin, and horizontal error bars repr esent the rms scatter for each bin of selected galaxy 
sample values. All measured values of absolute bias, ro, and 7 are given in Tables [T|6] 
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Fig. 13. — The projected correlation function, u}p(rp), for red (left) and blue (right) galaxy samples as a function of SFR threshold, 
limited to z<1.05. The dark matter projected correlation function is shown for 2=0.9, near the moan redshift of all galaxy samples. 



ing that they occupy the same mass halos at these red- 
shifts. The overall trends in clustering amplitude agree 
with z 1 environment studies performed in CP08, who 
found that the galaxy overdensity increases at both ex- 
tremes of galaxy SFR. 

Tables |3] and H] list the ro, 7, and bias values for each 
SFR sample, and Table [3] further lists the minimum and 
mean halo mass estimated for the SFR threshold sam- 
ples. We note that while there is a significant rise in the 
clustering amplitude on small scales for the highest SFR 
galaxies, there is no trend between SFR and 7, within the 
errors. The red galaxy SFR samples, which are just as 
clustered on large scales as the highest-SFR blue galaxy 
sample, have a consistent clustering slope as all of the 
blue SFR samples, though the errors on the slope are 
large. 

Interestingly, when com paring the stellar mass and 



SFR samples in Figure 12 the log(V')>1.5 Mq yr ^ blue 



galaxy sample has similar correlation lengths as the high- 
stellar- mass red ga laxy samples. Further, as shown in 
Figures [TT] and 14 we see similar enhanced clustering 
on smallscales relative to the best fit large-scale power 
law behavior in these galaxy samples. The similarities 
between high-SFR blue galaxies and high-stellar mass 
red galaxies on both large and small scales indicates that 
they are found in similar mass halos, perhaps linked to a 
shared star formation history, with high-SFR blue galax- 
ies evolving to high-stellar mass, low-SFR red galaxies 
at later epochs after the star formation in the galaxy is 
quenched. 

5.3. Clustering Dependence on Specific Star Formation 

Rate 

We also measure galaxy correlation properties as a 
function of specific SFR (sSFR). The sSFR is the in- 
verse of the star-formation timescale of a galaxy. In sec- 
ular evolution scenarios, the measured sSFR timescale is 
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Fig. 14. — The relative deviation in £jp(rp) from the best-fit power law for red (left) and blue (right) galaxy samples binned by SFR. 
Similar to stellar mass in Figure refirelmasscorrcolor, the power law is fit to u]p{rp) on scales of 1-10 Mpc, and each sample is offset 
by an arbitrary 0.3 dex in the figure for clarity. Wc find evidence of enhanced small-scale clustering in the highest-SFR blue galaxy sample 
but not for either red galaxy sample. 



proportional to the existing stellar mass relative to the 
amount of cold interstellar gas available for production 
of stars. Shorter timescales indicate rapid star formation 
from large, cold interstellar gas reserves and vice versa 
for longer timescales. As shown in Figure [2] sSFR is 
strongly correlated with restframe color as it reflects the 
star formation history of the aggregate stellar population 
within a galaxy. 

Figure 
threshold' 



15 



shows w„(rp) for the color-depe nde nt sSFR 



samples described in Section |3.3| The 
sSFR samples are complete between 0.74<z<1.05 for 
Mb<—20.5 blue galaxies and AfB<— 21.0 red galaxies. 
In general, lower sSFR samples are more clustered than 
higher sSFR samples, which is expected given the strong 
correlation between clustering and galaxy color (COS). 
Figure [T6| shows the deviation from the best fit power law 
for sarnples binned by sSFR. We find that for the lowest 
sSFR galaxies in both red and blue samples, ujp{rp) devi- 
ates significantly on small scales (rp<0.3 h^^ Mpc) rel- 
ative to the large-scale power law fit. There is a factor 
of 4 increase in the small-scale clustering amplitude for 
— 10.6<log(?/'/M*)<— 9.9 red galaxies and a similar in- 
crease for — 8.95<log(i/'/M*)<— 8.6 blue galaxies. The t- 
test for both these samples produce p- values of p = 0.013 
and p = 0.018, respectively. Additionally, we find that 
the 0.74<z<1.4 blue galaxy sample with A/b<— 21.0 and 
— 8.95<log(-0/Af*)<— 8.6 also deviates from the large- 
scale power law fit by an average factor of 7 and has a 
highly significant increase in clustering {p — 2.8 x 10^^) 
on rp<0.3 Mpc scales. The enhanced small-scale 
clustering indicates that the satellite fraction is related 
to how efficiently galaxies within a single halo, regard- 
less of restframe color, process gas and produce stars. We 
discuss the implications of this result further in Section[6j 
The right panel of Figure 12 shows the best fit tq vai- 
ues for galaxy samples binned by sSFR. The clustering 
length rises for blue galaxies with lower sSFR values and 
is constant for red galaxies with log(7/i/Af*)< — 10. Both 
red and blue galaxies with — lO.5<log(-0/M*)<— 9.5 have 



similar clustering lengths and therefore likely occupy sim- 
ilar mass dark matter halos. The correlation slope, 7, 
also roughly trends towards higher values for lower sS- 
FRs, although the trend is not highly significant given 
the errors. Tables [5] and |6] list the tq, 7, and bias val- 
ues for each sSFR sample. Table [5] further records the 
minimum and mean halo mass estimated for the sSFR 
threshold samples. 

6. DISCUSSION 

The primary goal of this work is to measure the cluster- 
ing properties of complete z ^ 1 galaxy samples selected 
by stellar mass, SFR, and sSFR . We present a macro- 
scopic view of our results in this section. To facilitate this 
discussion, we separate our findings into small- and large- 
scale behavior at rp ~1 Mpc, which roughly corre- 
sponds to the transition between the one and two-halo 
terms in HOD models. We first discuss the large-scale 
clustering results, including a comparison with previous 
work and relevant environment studies, and then discuss 
the small-scale clustering results. 

6.1. Large-scale Clustering Behavior 

The projected 2PCFs measured in this study allow us 
to quantify the relationship between galaxy clustering 
and stellar mass and SFR at 2; ~ 1. We first discuss how 
our power law fits to the projected correlation function 
on large scales compare with other studies. In the left 
hand panel of Figure |17[ we compare our ro results for 
color-independent , stellar mass thr eshold sampl es to re- 
sults f rom NMBS dWake et al |20Tl] ) and VVDS ( [Meneux 



et al. 



2008) at z ~ 1. 1 Wake et ai. calculated the angular 



ation function ot galaxy samples as a function of 
stellar mass using NMBS photometric redshifts in por- 



tions of the COSMOS and AEGIS fields iScoville et al. 



2007| [Davis et al.||2007[ [van Dokkum et al!pOn5| ), cover 
ing 0.4 deg^ in total. For NMBS, the angular correlation 
function is fit at 0.9<2;<1.3 over scales of ~ 0.05 — 10h~^ 
Mpc with a fixed 7 = 1.6. The DEEP2 results shown also 
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Fig. 15. — The projected correlation function, ajp(rp), for red (left) and blue (right) galaxy samples selected using sSFR thresholds 
between 0.74<2< 1.05. The dark matter projected correlation function is shown for z=0.9, near the mean redshift of all galaxy samples. 
The dotted line in the left panel corresponds to ajp(rp) for blue galaxies with log(i/i/M«)<— 8.6 yr~^, shown in the right panel. 
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Fig. 16. — The relative deviation in uip{rp) from the best-fit large-scale power law for red (left) and blue (right) galaxy samples binned 
by sSFR. A significant rise in the clustering amplitude is seen on small scales for the highest sSFR samples for both red and blue galaxies. 



have a fixed 7 = 1.6 to facilitate comparison and remove 



degeneracies. iMeneux et aL used spectroscopic r edshifts 
from the VVDS-Deep held ( |Le Fevre et al.||2003[ ) to cal- 
culate Wp(rp) over 0.49 deg'^, fitting vq over rp=0.1-21 
Mpc. While the VVDS results also allow 7 to vary 
in the fits, a slightly higher value of 7 = 1.8 is typical 
for their 2PCF data. We find good agreement with the 
VVDS clustering amplitudes; both VVDS and DEEP2 
results are less clustered than the NMBS results, though 
the difference is not significant given the NMBS errors. 
This may be due to cosmic variance within the NMBS 
z = 1.1 sample, where independent correlation functions 
in the COSMOS and AEGIS samples are quite different 
on large scales (although they statistically agree given 
the size of the errors, see their Figure 2). 

We also note that the DEEP2 stellar mass 
threshold sample with galaxies of all colors above 
log(M*/M0)>lO.5 is somewhat incomplete for red galax- 



ies at the low mass end beyond z>0.9. As stated 
in Section |3.1[ the red galaxy incompleteness between 
lO.4<log(M*7^/0)<lO.8 is estimated to be ~ 20% be- 
tween 0.74<z<1.05. As the clustering amplitude at a 
fixed stellar mass appears to be similar for red and blue 
galaxies above log(M*/M0)>lO.8, we expect the cluster- 
ing strength for the log(Af*/Af0)>lO.5 threshold sample 
to be relatively unaffected by the red galaxy incomplete- 
ness. If we assume that rg for the color-independent 
sample behaves as a weighted average between the red 
and blue galaxy samples (given in Table fl]), then in- 
creasing the number of red galaxies by 20%would only 
increase the clustering strength by ^ 1%, much less than 
the quoted error. 

We further compare our best fit clustering amplitudes 
to similar measurements in the local universe from SDSS 
(|Li et alj[2006) 



17 shows 



The center panel of Figure 
vq for color-independent, binned stellar mass samples in 
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SDS^ As |Li et al.| ( [2006[ ) did not measure tq from their 
projected 'AffJb's, we tit their ujp{rp) data over rp=l- 



10 h 

data. 



^ Mp c scales i n a sim ilar manner to our DEEP2 
2006) ujp{rp) data is drawn from 



The Li et al 

a sample of ~ 2U0, OCio bl3SS galaxies with an average 
statistical error of 5% per measured scale over the fitted 
scale range and does not account for the full covariance 
between measured scales. Therefore, our power law fits 
to their published data have a small 1% statistical error 
for each stellar mass sample and may not be fully re pre- 
sentative of the t rue error in the SDSS ujp^rp) data. Ze- 
havi et al. (2011) estimate the rg error to be ~5% from 
jackknited samples of galaxy luminosity bins in SDSS, 
and therefore our conclusions comparing clustering am- 
plitudes should be relatively robust. The same plot also 
shows the three most complete stellar mass samples in 
VVDS. Because the VVDS-Deep data becomes incom- 
plete for blue galaxies below log(MH,/M0)<lO, the mea- 
sured clustering amplitude i s systematically biase d to 
lower values of tq. At most, Meneux et al. (2008) cal- 
culated that the (9.5<log(M*/iWo)<10.0) stellax mass 
bin underestimates the clustering amplitude by 10% rel- 
ative to mock catalogs. DEEP2 color-independent stellar 
mass samples are limited to log(M*/M0)>lO.5, while the 
VVDS data extends to somewhat lower stellar mass with 
an /ab=24 magnitude limit. Again, we find excellent 
agreement in clustering strength between DEEP2 and 
VVDS at z ^ 1, both of which have significantly lower 
tq values than local SDSS galaxies for all stellar masses 
probed. At a fixed stellar mass, we estimate that the 
clustering amplitude in color-independent samples has 
increased by roughly 35% from z = 1 to z = 0.1. 

The right panel of Figure [17] shows rg for stellar mass 
selected samples separated by red and blue restframe col- 
ors in both SDSS and DEEP2. On large scales, we find 
that blue galaxies are somewhat more clustered (« 15%) 
locally when compared to z ~ 1, although the rg val- 
ues are completely consistent at the highest stellar mass, 
log(M*/M0)=11.25. Red galaxies at the same stellar 
mass, however, are much more clustered locally than at 
2^1. We interpret this to mean that lower stellar mass 
blue galaxies have much of their current halo mass in 
place by z ^ 1, and therefore there is little evolution in 
the clustering amplitude over the last 8 Gyrs. However, 
the halos that host massive red galaxies, which are more 
likely to be central galaxies, accumulate fractionally more 
halo mass relative to their existing stellar mass over the 
same span of cosmic time. This effectively causes the 
clustering amplitude to increase at a fixed stellar mass 
and therefore become more clustered on large scales at 
recent epochs. These results support stellar-halo mass 
assembly history models where lower stellar mass blue 
galaxies have most of their final halo mass in place at 
z ~ 1 but continue to add to their stellar mass, while 
higher stellar mass red galaxies form most of their final 
halo mass at z<l but have little evolution in their stel- 



lar mass (Zheng et al. 


2007 


Coupon et al.|20l2 


). 



(Jur clustering measurements also agree with the blue 

The SDSS and NMBS stellar masses have assumed a Kroupa 
IMF, while DEEP2 and VVDS stellar masses use a Chabrier IMF; 
the differences in IMF produce a negligible difference of 0.05 dex 
in stellar mass for our clustering comparisons. 



and red luminosity-dependent clustering results of COS, 
who found that red galaxies are more clustered than blue 
galaxies at z = 1 at the same luminosity (which corre- 
sponds to a lower stellar mass for blue galaxies compared 
to red) and less clustered than their local red counter- 
parts at the same luminosity (rg sa 5 Mpc at z ~ 1 
versus rg w 6 h^^ Mpc locally, see Figure 11 in C08). 
Further, the results in the lower red galaxy stellar mass 
bin are consistent with the results found in C08, given 
the errors. We find a slightly larger clustering difference 
between z 1 red galaxies select ed by stel lar mass and 
those at z ^ 0.1, using fits to the Li et al. Wp(rp) data. 
We note that our highest stellar mass red galaxy sample 
is a threshold sample limited by the total probed volume 
in DEEP2 volume. This upper stellar mass limit will 
not include the rarest, most massive galaxies that would 
be present in the larger SDSS volume within an equiv- 
alent stellar mass range. Another concern is that less 
massive star-forming galaxies that have been reddened 
by dust could be included in red galaxy samples, which 
may cause rg to be underestimat ed. Results from the 
PRIMUS survey ( |Zhu et al.|[2011[ ) show that such heav- 
ily reddened star-tormmg galaxies are rare at the bright 
end of the red sequence (< 10%) at intermediate redshift 
and therefore should have relatively little weight in the 
large-scale 2PCF amphtude. 

There are relatively few studies that have explored 
galaxy clustering at z ~ 1 as a function of stellar mass, 
SFR, and sSF R to which we can compare directly. How- 
ever, recently Lin et al. (2012) performed an angular 
clustering analysis of i^ziv -selected star-forming galax- 
ies at z ~ 2. They also found that the measured cluster- 
ing scale length increases strongly with increasing stellar 
mass (log(M*/M0)>9), increasing SFR (log(V')>0.5), 
and decreasing sSFR {\og{i{}/ M^)<—8.6), in good agree- 
ment with the trends found here. However, because 
we have restricted our blue galaxy sSFR samples to 
Mb<— 20.5 and log('0/M*)<— 8.6 for completeness to 
z<1.4, we do not probe higher values of sSFR beyond 
z>l in this study. We therefore cannot confirm the ex- 
istence of a turnover in the clu stering behavior at higher 
values of sSFR as presented in Lin et al. (2012). 

We can a lso compare our results as a function of sSFR 
to those of 



Li et al. 



( |2008[ ), who measured Wp(»'p) for a 
large sample ot bUbS galaxies with sSFRs ranging from 
-ll<l og('0/MH,) <— 9. To facilitate this comparison, we 



follow Li et al. and compute the best fit ^(r) power 
law from 10 jackknifed samples of our Wp(rp) data as a 
function of sSFR and normalize to a constant power law 
of ^(^)pow = h"^ Mpc)^-^'^. As we fit ujp{rp) over 
scales 1-10 Mpc in each jackknife sample, we evalu- 
ate the power law ^{r)fit for the unweighted mean values 
of rg and 7 at a fixed scale of 5 h^^ Mpc and plot the 



ratio of the 2 PCFs as a function of sSFR (see Figure 18 ) 
Compared to Li et al. ( 2008) (also aX Vp — 5 h~^ Mpc) 
we find good qualitative agreement with the trend seen 
for local clustering results. At log(V'/M*)>-9.5, there 
is a significant trend to lower clustering amplitudes at 
higher sSFR, while the large-scale clustering amplitude 
is relatively independent of sSFR at log(V'/A/*)<— 9.5. 
At log(?/'/Af*)<— 10.5 we do not see a rise in the large- 
scale clustering amplitude, as is seen locally. The trend 
of increased clustering at lower sSFRs should be expected 
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Fig. 17. — (left) Clustering scale length, ro, measured from color-independent stellar mass threshold samples in DEEP2 (black square), 
WDS data (green star) from Mcneux ct al. (2008|, and NMBS data (brown circle) from Wake ct al. (2011). Data point locations indicate 
the lower mass threshold limit imposed on each s'ample. The Wake et al. measurements use NIK photometric redshifts to compute the 
angular correlation function at 0.9<z<1.3. Both WDS and DEEP2 have a lower clustering amplitude than the NMBS at log(A/, )<11.0 
h~ Mq, although the errors from the 0.4 deg^ NMBS s urvey are relat ively large in comparison, (center) ro resulting from power law fits 
to color-independent stellar mass bin samples in SDSS ( Li et aL]|2006 I, WDS, and DEEP2. Fits are performed over rp=l-10 Mpc 
in DEEP2 and SDSS (purple diamond) while WDS is fit over scales of rp=0.1-21 Mpc. All points have similar mass bins sizes of 
0.4-0.5 dex, and the power law slope 7 is not fixed in the ujp{rp) fits. The DEEP2 clustering scale lengths agree well with those from 
WDS, both of which show lower clustering amplitudes at earlier times for all stellar masses probed, (right) Best fit ro values for restframe 
color-separated red and blue galaxy stellar mass bin samples in SDSS and DEEP2 for a fixed 7 = 1.6 over rp=l-10 Mpc. DEEP2 
mass bins range in size from 0.3 to 0.4 dex for blue and red galaxies, respectively. Blue galaxies in the local universe are somewhat more 
clustered at a given stellar mass than at ^ ~ 1, while massive red galaxies are more clustered locally. 
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Fig. 18. — The best fit power law for ^(r) measured from the 
DEEP2 sSFR samples on 1-10 h'^ Mpc scales (red circles). The 
real-space 2PCF is normalized to a power law of ^(r)pow = {r/5h~^ 
Mpc)~^'^ at a fixed scale of r = 5h~^ Mpc to match the aver- 
age scale for the DEEP2 data. |Li et alT] |2008 ) produced similar 
measurements in SDSS at equivalent scales (black crosses). The 
relative trend in sSFR is similar between 2: ~ 1 and local redshifts 
for log('!/i/M«)> — 10.5; lower clustering amplitudes are observed 
at progressively higher sSFR values. For log(i/i/Af«)<— 10.5, the 
2 ~ 1 galaxies are less clustered than their local counterparts at 
fixed sSFR. 



given the COS results, as sSFR is highly correlated with 
restframe color (see right panel of Figure [2]) and COS 
found an increased clustering amplitude witn redder col- 
ors. 

Here again, as with the highest stellar mass red sample 
(log(M*/M0)>ll.O), the sample with the lowest sSFR 
(log(V'/M*)<— 10.6) is a threshold sample which is lim- 
ited by the DEEP2 survey volume. As the clustering 



signal from the rarest galaxies cannot be measured from 
DEEP2 data, the clustering amplitude in this sample 
may be underestimated relative to larger survey volumes 
such as SDSS. The effect is expected to be minimal. 

6.2. Clustering Amplitude and the SFR-M^, 
Relationship 

While galaxy samples selected by SFR probe a wide 
range of stellar mass, stell ar mass and SF R are known 
to be correlated (e.g., |Noeske et al.||2007a ) such that the 
mea n ste llar mass increases with increasing SFR. In Sec- 
tion [5]T] we showed that the clustering amplitude tq in- 
creases monotonically with stellar mass, as the amount 
of baryons processed into stars correlates with the halo 
mass. Therefore, it is possible that the change in mean 
stellar mass for our blue SFR samples is responsible for 
the observed increase in clustering amplitude with in- 
creasing SFR^_SimilBi_jnv2S^ 

studies fPeng et al.]|2010[ [Sobral et a"L]|2010| ) have found 
little to no ditterence between the Sl^'R-density and stel- 
lar mass-density relationships, indicating that the SFR- 
density relation could be entirely due to differences in 
stellar mass. 

To address this question, we generate a prediction for 
ro based on the measured clustering with stellar ma ss us- 
ing a fixed slope of 7=1.6 (see left panel of Figure 12 ). We 
fit an exponential relation between tq and log{MZ/ Mq) 
for binned blue galaxy stellar mass samples, finding a 
best fit relation of 



„2.571og(M./M0)-27.1 



2.^ 



(9) 



We then calculate a prediction for tq as a function of 
SFR by weighting the fit relation with the stellar mass 
distribution of each SFR-selected sample. Because ^(r) is 
actually a pair-weighted statistic, this method is only an 
approximation to the tq that would be measured from a 
given stellar mass distribution. The ratio of the mea- 
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sured tq to the pre dicted rg as a function of SFR is 
shown in Figure 19 If the large-scale clustering am- 



plitude observed can be predicted given the stellar mass 
distribution of each SFR sample and the observed rela- 
tion between tq and stellar mass, ?'o('0)obs/^o(''/')fit would 
equal one. We find that the mean r-o('!/')obs/'"o ("0)151 over 
all SFR is 1.057±0.025; the deviation from unity has a 
p-value oi p — 0.076 and therefore is significant at the 
p < 0.1 level. We also find that there is a trend towards 
larger deviations from the mass-predicted tq at higher 
SFRs, suggesting that most, but not all, of the SFR - tq 
relationship can be explained by the relationship between 
SFR and stellar mass. 

While most of the clustering amplitude as a function 
of SFR may be explained by the tq-M* relation, it is 
interesting to investigate the SFR-dependent clustering 
behavior where the stellar mass distributions of SFR- 
selected samples are equal. In Figure [20j we plot SFR 
versus stellar mass for our parent DEEP2 sample with 
contours showing lines of constant log(?/'/M*). As we 
find that rg increases with decreasing sSFR (see Fig- 



ure [12), the clustering amplitude increases from the up- 
per left (high sSFR) to the lower right (low sSFR) in this 
space. To test whether high SFR galaxies are more clus- 
tered than low SFR galaxies with similar stellar masses, 
we construct a galaxy sample limited in stellar mass 
to 10<log(M,/Mo)<ll and sSFR to log(V'/M*)<-9.9, 
which is complete in the volume between 0.74<z<1.05. 
Thi s sa mple selection is shown with a dotted line in Fig- 
ure [20j We then perform a linear fit to the SFR- A/* de- 
pendence within this selected sample, weighting by the 
average errors of 0.3 dex in M* and 0.25 dex in SFR. 
The resulting fit (shown with a solid black line) roughly 
corresponds to the lo cation of the star-for ming "main 
sequence" ((MS) e.g., [Noeske et"ar||2007al ) for the se- 
lected sample. Separatmg the galaxy sample into popu- 
lations above and below the MS, we find that the mean 
stellar mass is log(M*/Mo)=10.36 for both populations, 
while the average SFRs are log(7/')=1.24 Mq yr~^ and 
log(V')=0.84 Mq yr~^, respectively. Comparing their 
clustering properties, we find galaxies above the MS have 
a clustering amplitude of ro=3.73±0.18, while galaxies 
below the MS have ro=4.36±0.21 and are therefore more 
clustered at a given stellar mass. This confirms that the 
tq-sSFR results found above are likely not driven by stel- 
lar mass differences between the samples. 

While the stellar mass and SFR values used in this 
study are not necessarily highly accurate on an individ- 
ual galaxy basis, they do represent the global average 
well. Any scatter in these values would wash out the 
observed SFR-M* relation and clustering correlations in 
this space. In particular, because the highest SFR galax- 
ies have the bluest restframe colors at a given IR lumi- 
nosity, it is possible that our color-based stellar masses 
underestimate the true stellar mass as measured in the 
IR (Weiner et al., in preparation). We have compared 
our color-M/L stellar masses to galaxies with measured 
iiT-band stellar masses both above and below the MS and 
find no statistically significant bias in the mass distribu- 
tions. 

The clustering dependence with SFR at a given stellar 
mass is particularly interesting in light of the differences 
in galaxy properties observed for galaxies above and be- 
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Fig. 19. — The measured clustering scale length (ro(i/')obs) fof 
blue galaxy SFR samples normalized to the scale length predicted 
by the correlation between ro and stellar mass. The predicted scale 
length (''"o(V')fit) is generated from an exponential fit to the blue 
galaxy stellar mass samples (Figure[l2]l and weighted by the stellar 
mass distribution of each SFR-selected sample. The weighted mean 
over all SFRs is ro(V')obs/^o('/')fit=t-057ib0.025, and a weighted 
linear fit (red, solid line) shows a trend for larger deviation from 
the stellar mass-predicted rg at higher SFR. 



low the MS. In general, galaxies above the MS are found 
to have higher star-forming surface densities, smaller 
sizes, more dust attenuation, and somewhat higher Sersic 
indices than galaxies on or below the MS ( Schiminovich 



et al. 2007 



Elbaz et al.||2011[ rWuyts et al.n20IT] F 

particular, both |Elbaz et ai.| ana |VVuyts et al.l interpret 
their results to support major mergers as a dominant 
mechanism of star formation quenching, where galaxies 
on the MS experience a merger event or some other in- 
stability, which leads to a central starburst and bulge 
formation, and eventually results in a quiescent ellipti- 
cal galaxy. In this picture, an individual galaxy would 
move from first being on the MS to being above the MS 
during the merger stage and then eventually move below 
the MS as the galaxy becomes quiescent (blue arrows in 
Figure [201). 



However, the clustering results presented here do not 
support this general picture, as galaxies above the MS 
are less clustered than those on or below the MS. There- 
fore the majority of the galaxy population above the 
MS (i.e. with a higher sSFR) can not have previously 
been on the MS, as the clustering amplitude can not 
decrease with time. Our clustering results instead sup- 
port a picture in which most galaxies initially lie above 
the MS, where they have higher SFRs at a given stel- 
lar mass due to internal (non-interacting) processes, and 
undergo secular evolution down onto and eventually be- 
low the MS as their star formation is quenched. While 
major merger events can certainly occur, their contri- 
bution to the overall clustering trends in field galaxies 
must be sub-dominant. The general evolutionary picture 
must move from the upper left in the SFR-stellar mass 
plane down to the lower right to match our clustering 
measurements. The clustering properties as a function 
of restframe optical color in both SDSS (Zehavi et al. 
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Fig. 20. — The dependence of SFR on stellar mass, overlaid with 
contours of constant sSFR. Dashed lines show the selection of star- 
forming galaxies with z<1.05 used to measure the clustering of 
galaxies above and below the star-forming "main sequence". The 
black solid line shows a linear fit to the SFR-stellar mass relation 
in this sample. Colored lines show two possible evolution scenarios 
for star-forming galaxies. The major merger scenario (blue arrows) 
predicts that galaxies evolve along the star-forming main sequence 
and undergo a merger event that initially boosts the SFR, followed 
quickly by a period of star formation quenching and progression 
onto the red sequence. The secular evolution scenario (red arrows) 
assumes that galaxies initially lie above the main sequence, during 
a period of intrinsically rapid star formation and are eventually 
quenched through internal or secular processes, moving through 
the main sequence in the process. We find that galaxies above 
the main sequence have a lower clustering amplitude than galaxies 
below the main sequence at a given stellar mass. As the clustering 
of individual galaxies increases with time, our measurements favor 
the secular evolution scenario, as the galaxy population above the 
main sequence can not be dominated by galaxies that previously 
had a lower SFR at a given stellar mass, which are more clustered. 



20111 and DEEP2 (COS) also support this scenario, as 
resttrame optical color is well correlated with sSFR at 
z ~ 0.1 and z ^ 1, and both studies find that within the 
star-forming population, bluer galaxies are less clustered 
than redder galaxies. However, the results presented here 
do not probe additional evolutionary scenarios whereby 
the large-scale environment or minor merger events may 
play a role in galaxy quenching. 

6.3. Comparison with Environment Studies 

An alternate way to study the dependence of galaxy 
properties on large-scale structure is through overdensity 
measurements of a galaxy sample relative to the mean 
density of galaxies at the same redshift. Such studies 
measure how a particular galaxy property is correlated 
with environment and are inherently limited to probe 
physical scales larger than the density-defining popula- 
tion. At 2^1, a limiting scale of > 1 Mpc is com- 
mon. Several environment studies at z~l have found 
a strong correlation between stellar mass and density, in 
that galaxies in higher dens ity environments t end to have 



al. 


2006 




2010 



ellar mass dKauffmann et al. 2004| iCucciati etj 
[Cooper et ai.pOO^l IPeng'eTaTP OlO: Sobral ef 



2010). Our results are in broad agreement with these 



studies, as we find higher clustering amplitudes within 
increasing stellar mass between 9.6<log(Af*/Af0)<ll.O, 
above which the clustering amplitude remains constant. 



There has been some disagreement in the literature, 
however, as to whether the color-density relationship 
evolves as a function of redshift. Several previous studies 
of galaxy environments at z '--^ 1 have failed to detect any 
statistically significant relationship between colo r and 



density, in contrast to studies at lower redshifts ( Cue 
ciati et al.| |2006l I Scodeggio et all [20091 |Grutzbauch"et 



^ On the other hand, Cooper et al.| ( [2010[ 

probed the most massive DEEP2 galaxies and found 
a statistically-significant correlation between color and 
galaxy overdensity, with high mass red galaxies found 
preferentially in the most dense environments at z ^ 1. 
[Cooper et al. emphasizes that while no significant trend 
in the color-density relation may be detected in several 
environment studies, it is possible that a color-density 
correlation may simply be missed due to larger statis- 
tical errors and/or observational systematics at higher 
redshifts (e.g due to low sample completeness in VVDS 
at z>l). In our analysis, we do not detect a difi^erence 
between the clustering of red and blue galaxies at an 
equivalent stellar mass. However, our samples are defined 
to probe the bulk color-density relationship at z ~ 1, 
and we do not probe the most over- or under-dense re- 
gions specifi cally, which was a primary focus of [Cooper 



ga 



et al. J 2010 ). Our results are consistent with Cooper et 
( 2010 ) m that the color-density trend for high-mass 
axies at z^l is weaker than the stellar mass-density 
trend; the typical densities of red and blue galaxies at a 
given stellar mass are not significantly different. 

Other environment studies have measured the SFR- 
density relationship as a function of redshift and stellar 
mass. Locally, this relation is monotonic, with lower-SFR 
galaxies f ound in higher density environments Gomez et 
al. (2003). At earlier times, the SFR-density relation- 



smp was quite different; dense enviro nments are associ 



ated with both low and h igh SFRs ([Elbaz et al.[[2007 
CP08, [Sobral et aLl[2010[ ). Our SFR-ciustering results 
at 2 ^ 1 are in broad agreement with these studies, 
with the highest-SFR blue galaxies having roughly the 
same clustering strength as the lowest-SFR red galaxies. 
The similar environments found for galaxy populations 
with vastly different SFRs (log(i/')< 0.5 Mq yr~^ and 
log(V^)>1.5 Mq yr^^) indicate that the processes leading 
to the large-scale clustering and quenching of star forma- 
tion in massive galaxies are primarily independent. One 
possible scenario is that rapidly star-forming blue galax- 
ies at these redshifts are proge nitors of the quiesce nt red 
galaxies found at later epochs ( Cooper et al.|2006 COS). 
The preponderance of quiescent galaxies found in high 
density regions at 2 ^ 1 and the relative paucity of high 
SFR blue galaxies at low redshifts gives the appearance 
that the large-scale environment and quenching are re- 
lated. 

Mechanisms proposed to explain the evolution of the 
observed SFR-density relation involve the quenching of 
star forma tion in a varie t y of isolated and merging sys- 
tems ( So bral et al.|[2010[ |Peng et al.[[2010[) and include 
effects such as virial heating ot cold gas needed for star 
formation and the rate of gas exhaustion. However, as 
pointed out in the previous section, it seems likely that 
quenching specifically from major mergers plays a sub- 
dominant role at 2 ^ 1; the clustering of star- forming 
galaxies at a fixed stellar mass increases monotonically 
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with decreasing SFR, not vice-versa. Our results sup- 
port the theory that a simple sSFR-M* relation, such as 
proposed "mass-quenching" models where the quench- 
ing rate for galax ies at or above M* is proportional t o 
their stellar mass ( |Noeske et aH2007b[[Pe"ng et al.|2010] , 
should reproduce the gross relationship between star- 
formation, stellar mass, and clustering on large scales. 
However, we caution that not all quenching in star- 
forming galaxies must follow this form, as our samples 
probe only the average clustering properties of the z ^ 1 
galaxy distribution. 

6.4. Small-scale Clustering Behavior 

We now turn to our clustering res ults on small scales 
below Tp ^1 h^^ Mpc. In Section 



^(rp,7r) for blue galaxies selected 



5.1 



we found that 
Dy stellar mass ex- 
hibits stronger "Fingers of God" on small scales at higher 
stellar mass (see Figure Isl), indicating that higher-stellar 
mass blue galaxies reside in halos with higher velocity 
dispersions and therefore likely reside in higher mass ha- 
los than lower-stellar mass blue galaxies. Such an obvi- 
ous increase in the "Fingers of God" within ^(rp, tt) is not 
seen for blue galaxy samples selected by SFR or sSFR, 
suggesting that these properties do not correlate with the 
dark matter halo mass as directly as stellar mass does. 
Although we do not show them here, we find that all 
red galaxy samples exhibit strong "Fingers of God" , in- 
dependent of the stellar mass or SFR and similar to the 
red luminosity samples in COS (see their Figure 8). 

In order to reveal any relative change in the small- 
scale clustering amplitude, we divide Wp(rp) on all scales 
by the best fit power law on scales 1-10 h^^ Mpc and 
look for excess clustering below rp<0.3 h^^ Mpc rela- 
tive to this best fit power law. In HOD models, the 
slope of a;p(rp) on these small scales, where the one- 
halo term dominates the correlation function, is governed 
by the relative number of central and satellite galaxy 
pairs within an individual halo, as well as the location 
of the probed halo mass scale on the halo mass func- 



tion (iB crlind fc Weinberg' 2002 ; ' ^havi et al.|2005{ [Zheng 
et al.T|2005, Skibba fc Shcth 200?| ). h'ov massive galaxy 



samples with a minimum halo mass corresponding to 
> 10^^'^ h~^MQ at z = 1, the Wp(rp) amplitude rapidly 
increases on small scales due to a greater contri bution 
from central-satellite a nd satelhte-sa tellite pairs (Sheth 
[fc To rmen (2002 ); Tinker et al.| ( [2008| ; see also Appendix 
A ot Zheng et aL '] 2009 ) tor further details) . Such mas- 
sive halos are near the exponential tail of the halo mass 
function, and higher-mass halos are increasingly rare in 
the volume surveyed by DEEP2. 

For our galaxy samples defined by stellar mass, 
we find no significant increase in clustering on small 
scales for any of the blue galaxy samples. How- 
ever, both mass-selected red galaxy samples exhibit 
enhanced clustering on scales below rp<0.3 Mpc 
{p = 0.005 for lO.5<log(M*/M0)<ll.O and p = 0.04 
for log(M*/A/0)>ll.O). This difference between blue 
and red galaxy mass samples may reflect a higher contri- 
bution from central-satellite pairs, due to both a lower 
satellite fraction for massive red galaxies (compared to 
blue galaxies) and the higher minimum halo mass (~ 
10^'^-^ h~^MQ), as discussed above. The lack of a de- 
tectable small-scale rise for the blue galaxy sample at an 



equivalent stellar mass may be due in part to larger mea- 
surement errors in t he highest mass sample (see the right 
panel of Figure 11) and/or a larger contribution from 



satellite-satellite pairs than central-satellite pairs, which 
smoothes the transition in cjp(rp) between the one-halo 
to two-halo terms. The small-scale clustering in color- 
independent stellar mass samples is driven primarily by 
the fraction of red to blue galaxies at a fixed stellar mass. 
Accordingly, we observe that the color-independent sam- 
ples exhibit a similar small-scale rise at high stellar mass 
due to the high fraction of red galaxies present in the 
sample. 

When binned by SFR, we find that red galaxies with 
low SFRs do not show a rise on small scales, while blue 
galaxies do exhibit a 7-fold increase in small-scale clus- 
tering amplitude over the best-fit large-scale power law 
below rp<0.3 h^^ Mpc. While the cjp(rp) error in the 
highest-SFR blue galaxy sample (log(-!/))>1.5 Mq yr~^) 
is the largest of all measured blue galaxy samples, the 
deviation from the large-scale power law has a p value of 
p — 0.038 and is statistically significant at the p < 0.05 
level. One possible explanation for the correlation be- 
tween high SFR and increased small-scale clustering am- 
plitude is that star formation is being triggered by galaxy 
interactions, eit her through merger events or tidal in- 



teractions (e.g. Barton et al. 20071. Late- type galax- 



ies in the blue cloud contain large gas reservoirs which 
can support enhanced levels of star formation triggered 
though close interactions of galaxy pairs. Enhanced IR 
luminosity, which is tightly correlated with star forma- 
tion, has been previously observed in extremely blue 
DEEP2 galaxies in close kine matic galaxy pair s {5rp = 50 
/i-i kpc) within th e EGS ([Lin et al.[[2007[). Further, 



Robaina et al. ( 2009 ) show enhanced clustering with SFR 
on slightly larger scales (40<r'p<180 kpc) between 



0.4 <z<0 .8 in a sample of COMBO-17 galaxies ( [Wolf et 
aL][2003 ) matched to Spitzer 24/im measurements. 

Our results show a similar increase in small-scale clus- 
tering power at high SFR, on scales rp<0.3 Mpc. 
However, we have demonstrated that an increased SFR 
is also correlated with higher stellar mass, and there- 
fore it is also possible that galaxies with higher stellar 
mass have a stronger one-halo term due to their higher 
halo mass. It is possible that both stellar mass and SFR 
enhancement are correlated with the increase clustering 
seen on small-scales in SFR-selected blue galaxies. We 
note, however, that red galaxies with low SFR have simi- 
lar halo masses to the high-SFR blue galaxies and do not 
show a similar rise, indicating that SFR enhancement is 
likely more strongly correlated with the increased clus- 
tering signal. 

These small-scale clustering trends allow us to inter- 
pret the behavior seen as a function of luminosity in COS 
(their Figure 10), who detected a relative small-scale rise 
for the brightest blue galaxies and a marginal increase 
for bright red galaxies. Figure [2] shows that blue galaxies 
with the highest SFRs also have the highest luminosities, 
and therefore the observed small-scale rise for bright blue 
galaxies may have been due in part to both higher stel- 
lar masses and SFR enhancement. The small-scale rise 
observed here for massive red galaxies is likely reflected 
in the marginally-significant increase seen for bright red 
galaxies in COS. As the mass-to-light ratio varies across 
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the red sequence, a given luminosity range corresponds 
to a wider range in stellar mass; if stellar mass is more 
strongly linked to halo mass, luminosity-selected samples 
will have a broader host mass range and hence exhibit a 
weaker small-scale clustering signal. 

Turning to sSFR, we find that galaxies that are form- 
ing stars most rapidly relative to their existing stellar 
mass have an increased clustering amplitude on small 
scales. Interestingly, both the high-sSFR sub-samples of 
the blue and red galaxy population show a small-scale 
rise at rp<0.3 Mpc. In each sample, we measure a 
4- fold increase in ojp{rp) relative to the large-scale power 
law fit over rp=l-10 Mpc scales. The increase in 
clustering is significant at the p < 0.05 level in both 
galaxy samples. Unlike the previous small-scale cluster- 
ing results with SFR, the sSFR samples have normalized 
the SFR relative to the stellar mass, and therefore the 
measured small-scale rise will be more directly connected 
to enhanced star formation at a fixed stellar mass. Be- 
cause sSFR is highly correlated with star formation his- 
tory and restframe color, the galaxies demonstrating this 
increased small-scale clustering are the bluest members of 
their respective populations (e.g. the bluest blue galaxies 
a nd the bluest re d galaxies). 



Li et al. (2008) measured the clustering of SDSS galax- 
ies selected by sSFR and found enhanced small-scale 
galaxy clustering in the local Universe, with >40% of 
these galaxies having close {vp <100 kpc) companions. 
They conclude that the observed clustering behavior is 
a signature of tidal interactions between galaxy pairs in- 
side the same dark matter halo, which leads to an in- 
fiow of cold gas and enhanced star formation. Evidence 
for enhanced sSFR on small scales was also measured 
from close pairs of blue galaxies (rp<50 h^^ kgc) at 



higher reds hift in PRIMUS between 0.25<2<0.5 (|Wong 
et aI]|20 lT1 ) and in DEEP2 betwe en 0.1<z<l.l 



;in et 



ar' 2D07p . Robaina et al. (20091 also found enhanced 
clustermg on scales rp<40 kpc tor star- forming galaxies 
with M* > 10^° Mq between 0.4<z<0.8. While we de- 
fine "small-scales" at a larger scale threshold of rp<0.3 
Mpc, we also find statistically significant evidence 
for enhanced small-scale clustering with sSFR. These re- 
sults support our interpretation that enhanced star for- 
mation, possibly triggered from galaxy tidal interactions, 
is correlated with the observed small-scale rise in our blue 
galaxy sample with high sSFR. 

However, enhanced small-scale clustering for high- 
sSFR red galaxies may be somewhat surprising. The 
range of sSFR for which we find this enhancement 
(— lO<log(-0/M*)<— 11) is often referred to as the "green 
valley" , a transition region fr om blue star-forming galax- 



ies to red quiescent galaxies ( 


Martin et al.|2007 


Salim et 


al. 


2009 


Mendez et al. 


2011 


). I'he enhanced c 


lustering 



could be due to recent merger events or tidal disruptions 
between red and blue galaxy pairs where the blue galaxy 
has provided ad ditional gas for star formation in the qui- 
escent galaxy. |Robaina et al. (2009) found that major in- 
teractions can contribute significantly to the small-scale 
clustering amplitude for the most dust-obscured star- 
bursts, which may be selected in our DEEP2 sample 
as red with the restframe color cut of Equation [T] Al- 
ternatively, close interactions between blue galaxy pairs 
could temporarily boost the SFR and eventually lead 



to quenching, which would move galaxies from the blue 
cloud to the red sequence. As we argued in the last sub- 
section, however, such events must be sub-dominant to 
the general quenching trend in individual halos. 

7. CONCLUSIONS 

In this study, we measure the two-point correlation 
function of complete DEEP2 galaxy samples selected by 
stellar mass, SFR and sSFR and separated by restframe 
color. We fit a power law to the projected correlation 
function and measure the bias and clustering scale length 
and slope on scales of 1-10 Mpc (fp=4.1 h^^ Mpc), 
where the two-halo term dominates the 2PCF. We also 
study the relative shape of the correlation function on 
small and large scales, which depends on the satellite 
fraction and dark matter halo mass. We summarize our 
findings as follows: 

1. The large-scale galaxy clustering amplitude 
increases monotonically as a function of 
stellar mass for star-forming galaxies with 
9.6<log(M*/M0)<ll, indicating that stellar 
mass closely tracks halo mass at z ~ 1. Within 
the limited mass range probed for red galaxies 
(lO.5<log(Af*/M0)<11.5), we find no signifi- 
cant trend between stellar mass and clustering 
amplitude. Within the measurement errors, our 
results agree with a z ^ 1 clu stering analysis from 



VVDS ( Meneux et al. 2008[ ) in the stellar mass 
range where both UEEP2 and VVDS samples are 
complete. 

2. Within the blue, star-forming galaxy population, 
the large-scale clustering amplitude increases as a 
function of increasing SFR. Red galaxies with low 
SFRs are strongly clustered and show no depen- 
dence on SFR within the range probed. We find 
the highest-SFR blue galaxies have the same clus- 
tering amplitude as red galaxies, as seen in previous 
DEEP2 environment studies. 

3. The observed trend between large-scale clustering 
amplitude and SFR can primarily be accounted for 
by the observed correlation between clustering am- 
plitude and stellar mass. However, the ro-stellar 
mass relation alone does not fully predict the clus- 
tering amplitude at all SFRs; there is a small clus- 
tering excess at high SFR above what is predicted 
by the stellar mass of each sample alone. This sug- 
gests that most, but not all, of the correlation be- 
tween large-scale clustering amplitude and SFR can 
be attributed to the SFR-stellar mass relation. 

4. Galaxy samples selected by sSFR show similar 
large-scale clustering trends as samples selected 
by stellar mass. The clustering amplitude de- 
creases with increasing sSFR (corresponding to 
bluer restframe color) while red galaxies with 
low sSFRs have a constant clustering strength in 
the range probed. While the large-scale cluster- 
ing amplitude only mildly depends on sSFR for 
— 11.0<log('(/'/-M*)<— 9.5, To drops significantly be- 
tween — 9.5<log(V'/M*)<— 8.5, similar to the trend 
observed locally in SDSS. 
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5. By constructing stellar mass-limited star-forming 
galaxy samples above and below the SFR-stellar 
mass "main sequence" , we find that the clustering 
amplitude increases with decreasing SFR at a fixed 
stellar mass (i.e. below the main sequence), which 
confirms that the clustering trends seen with sSFR 
are not driven entirely by stellar mass. Given that 
galaxies below the main sequence arc more clus- 
tered than those above, the bulk of the population 
above the main sequence can not be dominated by 
galaxies that used to be on the main sequence and 
are currently undergoing major merger events. In- 
stead, galaxies must smoothly evolve from above 
the main sequence; to below it, as the clustering of 
galaxies can only increase with time. 

6. We detect enhanced clustering on scales less than 
rp<0.3 Mpc, relative to larger scales, for high 
stellar mass red galaxies, high SFR blue galaxies, 
and the highest sSFR sub-samples of both red and 
blue galaxy populations. The increased small-scale 
clustering may reflect a combination of effects, in- 
cluding a changing satellite fraction, higher halo 
mass, and/or enhanced SFR. Wc conclude that 
triggered star formation due to close galaxy in- 
teractions is a likely explanation for the enhanced 
clustering seen at high sSFR. 

As our analysis uses the two-point correlation func- 
tion formalism and measures clustering properties 
for several important galaxy properties, we expect 
these results to be highly useful for galaxy evolution 
models. In particular, HOD and abundance matching 



models should greatly benefit from measurements 
of the clustering amplitude at z 1 with respect 
to stellar mass, SFR, and sSFR. We also anticipate 
that the bias values calculated here as a function of 
stellar mass and SFR will be used in projections of 
future BAO surveys that will probe the z>l universe. 
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TABLE 2 

Results for Stellar Mass Bin Samples 
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TABLE 3 

Results for Star Formation Rate Threshold Samples 
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TABLE 4 

Results for Star Formation Rate Binned Samples 
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1.01 


3.91 ± 0.27 


1.62 ± 0.15 


1.48 ±0.07 


Blue 


1.25-1.50 


1.35 


3243 


5.2 ±0.2 


3.2 


0.74-1.40 


1.11 


3.93 ±0.18 


1.48 ±0.22 


1.63 ±0.07 


Blue 


>1.50 


1.62 


1114 


2.3 ±0.1 


0.4 


0.74-1.40 


1.13 


4.45 ± 1.19 


1.38 ±0.27 


1.97 ±0.32 


Red 


-0.1-0.4 


0.14 


1081 


2.9 ±0.6 


2.6 


0.74-1.05 


0.87 


5.00 ± 0.29 


2.10 ±0.40 


1.72 ± 0.08 


Red 


>0.4 


0.65 


954 


4.7 ±0.3 


1.7 


0.74-1.05 


0.90 


5.26 ±0.63 


1.65 ±0.31 


1.90 ±0.13 



TABLE 5 

Results for Specific Star Formation Rate Threshold Samples 



Color 


sSFR 


<sSFR> 


^gal 


n 


Ccv 


z Range 


<2> 


ro 


7 


Bias 




<iW"haIo> 




[log 


:(yr-i)] 




[10-4 Mp 


e-3] 






[/i-i Mpc] 






[log(h 




Blue" 


<-8.60 


-9.13 


5777 


28.9 ±0.7 


8.5 


0.74-1.05 


0.91 


3.97 ±0.16 


1.47 ±0.08 


1.52 ±0.06 


12.07 


12.64 


Blue 


<-8.95 


-9.30 


3819 


19.2 ±0.5 


8.2 


0.74-1.05 


0.91 


4.43 ± 0.20 


1.63 ±0.07 


1.59 ±0.07 


12.20 


12.74 


Blue 


<-9.30 


-9.54 


1651 


8.2 ±0.4 


4.7 


0.74-1.05 


0.91 


4.66 ± 0.35 


1.66 ±0.14 


1.67 ± 0.09 


12.34 


12.84 


Blue' 


<-8.60 


-8.79 


7144 


13.1 ±0.4 


3.2 


0.74-1.40 


1.11 


4.21 ±0.11 


1.70 ±0.07 


1.64 ±0.05 


12.29 


12.80 


Blue 


<-8.95 


-9.11 


4693 


9.3 ±0.3 


2.9 


0.74-1.40 


1.07 


4.54 ±0.10 


1.82 ±0.07 


1.70 ± 0.04 


12.38 


12.88 


Blue 


<-9.30 


-9.54 


2013 


3.8 ±0.2 


1.5 


0.74-1.40 


1.03 


4.62 ± 0.24 


1.84 ±0.18 


1.76 ±0.08 


12.47 


12.94 


Red'' 


<-9.90 


-10.59 


1641 


8.2 ±0.4 


3.8 


0.74-1.05 


0.89 


5.17 ±0.23 


1.80 ±0.13 


1.84 ±0.06 


12.58 


13.02 


Red 


<-10.6 


-10.91 


795 


4.1 ±0.3 


2.6 


0.74-1.05 


0.89 


5.04 ± 0.32 


2.04 ± 0.25 


1.81 ±0.11 


12.54 


12.99 



" 2<1.05 galaxy samples are limited to Ma<-20.5. 
2<1.4 blue galaxy samples are limited to Ms<-21.0. 



TABLE 6 

Results for Specific Star Formation Rate Binned Samples 



Color 


sSFR 


<sSFR> 




n 




z Range 


<z> 


ro 


7 


Bias 




[log(yr- 






[10-4 : 


IVIpc-3] 






[/i-l Mpc] 






Blue" 


-8.60 - -8.95 


-8.80 


1959 


9.6 ±0.4 


3.2 


0.74-1.05 


0.91 


3.64 ± 0.36 


1.59 ±0.09 


1.42 ±0.11 


Blue 


-8.95 - -9.30 


-9.11 


2170 


11.1 ± 0.4 


3.8 


0.74-1.05 


0.91 


4.12 ± 0.43 


1.54 ±0.16 


1.54 ±0.10 


Blue 


<-9.30 


-9.54 


1651 


8.2 ±0.4 


5.0 


0.74-1.05 


0.91 


4.70 ±0.35 


1.62 ±0.13 


1.67 ±0.09 


Blue*- 


-8.60 - -8.95 


-8.79 


2451 


3.9 ±0.2 


2.9 


0.74-1.40 


1.12 


3.29 ±0.35 


1.30 ±0.28 


1.59 ±0.15 


Blue 


-8.95 - -9.30 


-9.11 


2681 


5.2 ±0.2 


1.4 


0.74-1.40 


1.07 


4.13 ± 0.21 


1.66 ±0.16 


1.64 ±0.07 


Blue 


<-9.30 


-9.54 


2013 


3.9 ±0.2 


1.6 


0.74-1.40 


1.03 


4.64 ±0.26 


1.83 ±0.17 


1.76 ±0.08 


Red" 


-9.90 - -10.6 


-10.28 


847 


4.2 ±0.3 


2.0 


0.74-1.05 


0.90 


4.82 ± 0.35 


1.43 ±0.32 


1.86 ±0.13 


Red 


<-10.6 


-10.91 


795 


4.0 ±0.3 


2.7 


0.74-1.05 


0.89 


5.04 ±0.32 


2.04 ±0.25 


1.81 ±0.11 



z<1.05 galaxy samples are limited to Mb<-20.5. 
2<1.4 galaxy samples are limited to Mb<-21.0 . 



